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I. 



INTRODUCTION 



Operations researchers, military planners and programmers, 
statisticians, economists, marketing personnel, managers and 
others are often faced with the need to analyze data in the 
form of a time series, which can be thought of as a sequence 
of observations of either a deterministic or stochastic 
process. In most cases, the objective of the analysis is 
to determine patterns or other recognizable behavior appar- 
ent in the data, and to then formulate a suitable mathemati- 
cal model for the time series from which forecasts of future 
behavior can be obtained. The value to a decision maker 
of being able to predict the future with some reasonable 
and statistically quantifiable degree of accuracy can not 
be overstated. For example, in areas such as budget expen- 
ditures, recruiting performance, commodity prices, population 
levels, resource consumption, manpower levels and consumer 
demand for products, decision makers charged with planning 
for the future should base their decisions at, least in 
part on the best available predictions about the future 
behavior of the time series in question. 

Until the late 1960 's, the techniques employed in time 
series analysis were primarily those of spectral analysis, 
with applications of harmonic analysis and mathematical 
transform theory. For an introduction to these topics, the 
reader is referred to Anderson [Ref. 2], Lewis [Ref. 7] 
and Jenkins [Ref. 10]. Due to the high degree of mathematical 
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sophistication required by the spectral analysis approach, 
the capability to perform time series analysis resided nearly 
exclusively with mathematieians and electrical engineers. 

As a consequence, the majority of decision makers came to 
use more understandable (and far less powerful) methods 
such as simple moving averages or exponential smoothing. 
However, in the late 1960 's and early 1970' s, the statistical 
analysis of time series by the methods developed by Box 
and Jenkins [Ref. 4] has gained widespread acceptance. 
Although these methods are mathematically nearly equivalent 
to the spectral approach (there are transformations that 
interconnect the two methodologies) , the Box- Jenkins approach 
is described in vocabulary more familiar to operations 
researchers, statisticians, economists and managers, and 
is therefore being used more and more by them in building 
models from measurements of the past. 

Many algorithms and computer programs for performing 
the analyses required by the Box-Jenkins approach have been 
developed and are available from several sources. Perhaps 
the best source is the collection of FORTRAN computer sub- 
routines which resides in the International Mathematical 
and Statistical Library (IMSL) [Ref. 9]. The primary prob- 
lem with using the available computer resources lies not 
in any deficiency of the programs themselves, but with the 
very nature of the Box-Jenkins methodology itself. The 
Box-Jenkins method is an iterative approach to time series 
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analysis, which is briefly described in Figure 1. [See 
Wheelwright and Makridakis, Ref. 16.] 

As indicated by Figure 1, the Box- Jenkins method is a 
multi-stage, iterative process. It begins with the postu- 
lation of a general class of time series models which has 
been found, experimentally, to be extremely rich. There- 
after, the procedure continues as a trial-and-error process, 
with decisions points where the analyst is required to 
select the next direction based on the best information 
available to him. Since each stage of the process outlined 
in Figure 1 may consist of several sub-steps, the modeling 
process itself can become quite time-consuming, even with 
the ready availability of the IMSL software resources. For 
example, an analyst employing the IMSL subroutines in a 
batch processing computer system to perform Box-Jenkins 
modeling of a time series might perform the following 
sequence of tasks : 

1. Prepare the time series data in the proper format. 

2. Plot and visually examine the time series, checking 
for nonstationarity , trends, seasonality, patterns, 
etc . 

3. Write a program to call the IMSL subroutine that 
will calcuate the mean, variance, autocorrelations 
and partial autocorrelations of the series. 

4. Plot the autocorrelations and partial autocorrela- 
tions; this provides much of the information re- 
quired for identifying the correct model for the 
series. 
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FIGURE 1. The Box- Jenkins Forecasting Methodology 
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5 . Write a program to call the IMSL subroutine which 
transforms the time series to adjust for seasonal 
patterns, nonstationary behavior or other behavior 
which deviates from that assumed for the class of 
models postulated. 

6. Repeat steps (2) through (4) with the transformed 
data. 

7. Review the statistical properties of the auto- 
correlations and partial autocorrelations for 
tentative model identification. 

8. Write a program to call the IMSL program that 
estimates the model parameters and computes the 
residuals . 

9. Write a program to perform goodness-of-f it tests 
for the model. 

10. Analyze the model residuals, using steps (1) through 
(9), as with the original time series. 

11. Refine the model using information gained by 
examination for any structure in the residuals. 

12. Repeat the preceding process until no structure 
remains in the residuals. 

13. When an adequate model has been obtained, write a 
program to call the IMSL subroutine which forecasts 
and determines confidence intervals for future 
values of the time series . 

Between each successive pair of steps, the analyst must 
manually intervene in the analysis process, making a decision 
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subjectively based on the information then available. 

Thus, much analyst interaction is necessary in order to 
determine a suitable mathematical model and forecast equa- 
tion. Even with rapid computer job turnaround time, the 
process described above could easily consume a full working 
day or more. 

This paper describes an effort to alleviate some of the 
organizational problems inherent in the modeling of time 
series using the Box-Jenkins iterative approach. This 
effort consists of an interactive computer program package 
which provides access in a structured way to the most useful 
IMSL and other subroutines necessary for Box-Jenkins time 
series analysis. This package, called the Time Series Editor, 
is written for use with the Naval Postgraduate School's 
Control Progam/Cambridge Monitor System (CP/CMS) . Since all 
working programs with the exception of the executive routine 
itself are written in FORTRAN, the Time Series Editor should 
be adaptable for use on other FORTRAN-capable time-sharing 
systems. The Time Series Editor assists the analyst in 
data preparation and entry, model construction and diag- 
nostic testing, and time series forecasting. In fact, 
with the User's Guide provided as an Appendix to this report, 
a complete Box-Jenkins time series analysis can be per- 
formed in a short time by even a novice computer user. 

For this reason, the Time Series Editor could be a valuable 
instructional tool for laboratory use in a time series 
course. 
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A brief overivew of Box-Jenkins methodology is provided 
in Chapter II, to serve as a point of reference for the 
material which follows. Chapter III contains a description 
of the algorithm employed to estimate non-linear least- 
squares parameters for a generalized (seasonal, nonseasonal, 
stationary, or nonstationary) Box-Jenkins model. Chapter IV 
contains descriptions of each of the major programs in 
the Time Series Editor. Chapters V and VI contain examples 
of actual use of the Time Series Editor, using a non- 
seasonal and a seasonal series, respectively. Chapter VII 
includes a summary of the report, and recommendations for 
future additions to the Time Series Editor. The Appendices 
include a User's Guide to the Time Series Editor, with 
instructions, sample user sessions and sample outputs, as 
well as listings of the Editor itself and all primary 
subprogram packages . 
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II. BOX- JENKINS METHODOLOGY: A PRIMER 



This chapter provides a brief overview of the methodology 
developed for the analysis of time series data by G.E.P. 

Box and G.M. Jenkins, in order to facilitate understanding 
of the models and programs employed in the Time Series 
Editor. For more detail on the material covered in this 
chapter, the interested reader is referred to the texts by 
Anderson [Ref. 1], Box and Jenkins [Ref. 4], Pindyck and 
Rubinfeld [Ref. 15] , Nelson [Ref. 14] , and Mabert [Ref. 12] . 



A. LISTING OF NOTATION 

The following is a listing of the notation employed 
in this chapter and in the remainder of this report. Each 
symbol is provided with a brief title or explanation of its 
use. 



U,} 

N 

- -2 
z, a„ 






{w^} 



a time series 

the length of the time series (number of terms) 

estimates of the mean and variance of the time 
series {z^} 

I ^ 

transform of {z^}, where z^' = (z^ + C) 
or z^ = ln(z^ + C) 

differenced series from {z }, where 

W = V^V^z' 
t St 

autocorrelation of lag k 

estimate of autocorrelation at lag k 

backward difference operator of order d 
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S 


seasonal backward difference operator of order 
D with season s 




backward shift operator of order s 


s 


length of the time series' season 




non-seasonal^autoregressive parameters 
(estimates, $) 


$ 


seasonal autoregressive parameters (estimates, $) 


9 


non-seasonal^moving average parameters 
(estimates, 9) 


0 


.A. 

seasonal moving average parameters (estimates, 0) 




series of white noise terms of the process 


2 

^a 


variance of the white noise process 


S(J,6) 


1 

sum of squared residuals of the model 


P 


nxamber of non-seasonal autoregressive parameters 


P 


nuiT±)er of seasonal autoregressive parameters 


q 


number of non-seasonal moving average parameters 


Q 


nxamber of seasonal moving average parameters 


’i^(B) 


transfer function of a linear filter 


'^kk 


partial autocorrelation of log k 



B. TIME SERIES, STOCHASTIC MODELS AND STATIONARITY 

A distinguishing feature of modern time series analysis 
is that the sequence of observations of a given variable is 
considered to be a realization of jointly distributed 
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random variables, such that the time series itself can be 
viewed as a stochastic process. A discrete stochastic 



time series, then, is a set of observations (z^ , Z2 , • • . , z^) 
generated sequentially in time by a set of jointly distributed 
random variables, where the actual data (z^, Z2, • . . ,z^) 
represents a particular realization of some joint probability 
distribution function f (z^^, Z2, • • • , z^) . With this function 
determined in some way, a forecast can be thought of as 

having been generated by a conditional probability distri- 
bution f (z^_|_j^ I Zj^, Z2, • . . / z^) . Now, the stochastic process 
which generates the time series (z^ , Z2 ^ ♦ f z^) is said to 
be stationary if its properties are unaffected by a change 
in time origin; that is, the process has reached a particular 
state of statistical equilibrium. This implies that the 
joint probability distribution associated with the m 

observations (z ,z ,...,z ), made at any set of times 

^1 ^2 ^m 

(tf, tz/ • • • / tm) » is identical to that distribution associated 

with the m observations (z^ +k'^t +k'‘‘"'^t +k^ made at 

^1 ^2 m 

times (t, +k, t_+k, . . . , t +k) . 

12 m 

For stationary discrete time series, the marginal 
probability distribution p(z^) is the same for all times 
t. Hence, the stochastic process embodied in the time 
series has a constant mean 



00 

y = E[z^] = / zP(z)dz , 

— CO 
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which defines the level about which it fluctuates in time, 
and a constant variance 



= E[ (z^ - U) ] = 



/ (z - ) p(z)dz 



which measures its variability about the level of the mean. 
Since the probability distribution p(z) is the same for all 
times t, its form can be discovered by plotting a histogram 
or relative frequency plot of the set of observations 
^^1' ^2' ’ * * ' ^t^ * Additionally, the mean and variance of 
the stochastic process can be estimated by the sample 
averages taken over time : 



N 




t=l 



and 



-2 

a 

z 



N 



- I 

N ^ 



(z^ ■ 



z) ^ 



t=l 



C. THE BASICS OF DIFFERENCE OPERATORS AND LINEAR FILTER 
MODELS 

This section provides a brief discussion of the differ- 
ence operators currently employed in writing Box- Jenkins 
models. Those most commonly used are described below in 
Table I. 
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Operator 



Name 



Definition 



B 



m 



Backward Shift Operator 
of order m 




z 



t-1 



z 



t-s 






V 



D 

s 



Forward Shift Operator 
of order m 



Backward Difference 
Operator of order D 
and season s 



F = b“^ 
Fz . = z 



t+1 



= 24 . , 
t t+m 



7'^ = (1 - E)° 



7“ = (1 - 



- ^t-1 



Table I. Basic difference operators. 



It is noted here that ^ 2 ^ could be written as a 



binomial expansion; that is. 



V°z^ = (1 - B)°z^ 



= ( I (°) (1)^ (-B)°"^)z^ 



j=0 



I (°)(-l)°-^z^_ 



D+j 



j = 0 



The stochastic models employed in the Box-Jenkins method 
are based on the idea that a time series in which successive 
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values are highly dependent can be regarded as having been 

generated from a series of independent "shocks” a^. These 

shock terms, sometimes referred to as white noise terms, 

are random drawings from some fixed distribution, generally 

2 

assumed Normal with mean zero and variance a . The white 

a 

noise process a^ is supposed transformed to the process z^ 
by what is known as a linear filter, shown in Figure 2 below 
[Ref. 4, p. 8]. The linear filtering operation takes a 



White 

Noise 






Linear 




Filter 





Process 
z . 



Figure 2. Representation of a time series as 
the output from a linear filter. 



weighted sum of previous observations, such that z^ can be 
written 



2^ = y + a^_3_ + a^_2 + 



= y + a^ . 



In general, y is the parameter that determines the "level" 
of the process, and 



21 



4^(B) 



2 

= 1 + * • * 

is the linear operator that transforms a^ into z^, called 
the transfer function of the filter. The sequence 

... formed by the weights may be (theoretically) 
finite or infinite; if this sequence is finite, or infinite 
and converges, the filter is said to be stable and the 
process stationary. The parameter y is then the value 
about which the process varies. Otherwise, is called 
nonstationary, and y has no specific meaning except as a 
reference point for the level of' the process. 

D. HOMOGENEOUS NON -STATIONARY AND SEASONAL PROCESSES 

In many cases of interest to the time series analyst, 
the process being examined is not stationary. Instead, 
the probabilistic structure of the process generating the 
time series changes with time. For example, there may be 
a trend or seasonal pattern in the time series. However, if 
the series can be found to exhibit behavior which is somewhat 
homogeneous, then the series can frequently be transformed 
into a stationary series that can be described by Box-Jenkins 
models . 

A series is said to be homogeneous nonstationary of 
order d if w^ = V^z^ is a stationary series for some integer 
d. Here V = (1 - B) denotes the difference operator such 
that V^a^ = The fact that a given series is 

nonstationary can be determined quickly by examination of 
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the autocorrelation function, to be discussed later. The 
series to be modeled should be differenced until the re- 
sulting series appears stationary or until the procedure 
appears to be making no improvement. 

A time series is called seasonal when it exhibits 
cyclical behavior over time, showing a regular period. 

Examples of seasonal processes include monthly rainfall, 
monthly crop yields, sale volume of bathing suits, livestock 
production rates and energy consumption. Seasonal patterns 
are often easy to spot simply by observing a plot of the 
series, or through knowledge of the process that generated 
the series. However, in many cases where the variability of 
the series is large, seasonal patterns will be difficult to 
distinguish from other fluctuations. Recognition of season- 
ality in a series is important, since it provides information 
useful in modeling and forecasting; as will be shown, the 
autocorrelation function can assist in recognizing seasonality. 
The Box-Jenkins modeling approach for seasonal nonstationary 
time series of season length s is to first transform the 
series using an appropriate seasonal differencing operator 
to "sweep out" the seasonal effect; that is, take 



If there is a trend or other types of non-stationarity pre 
sent, as well as the seasonal effect, the series can be 
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differenced again with the ordinary difference operator 
\intil stationarity is achieved. This transform would be 
written 



„ _ „d „D 

«t - V 2^ 



The transformed series is then modeled as a stationary 
series. For seasonal series where complete stationarity 
can not be achieved, a class of seasonal models is avail- 
able; these are discussed later in this Chapter. 

E. GENERAL CLASSES OF BOX- JENKINS MODELS 

This section will provide brief descriptions of several 
of the most commonly used classes of Box-Jenkins models, 
including autoregressive (AR) models, moving average (MA) 
models, mixed and integrated autoregressive-moving average 
(ARMA/ARIMA) models, and a generalized seasonal ARIMA model. 

1. Autoregressive Models 

In autoregressive (AR) models , the current value of 
the series under study is expressed as a linear combination 
of previous series values that explain the current obser- 
vation, plus an unexplained random (white noise) term a^. 

For a stationary series, or one that has been transformed 
to stationarity, the deviation from the mean value for each 
period, ~ modeled as dependent on weighted 

values of the previous deviations from the mean. The AR 
model can be written 
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^ *p^t- 



+ a. 



p t-p t 



In the expressions which follow the tilde will be dropped, 

but it is to be understood that the mean has already been 

subtracted from each observation. This equation is called 

an autoregressive process of order p, where (j)^ represents 

ttl 

a scalar weighting coefficient for the j previous period. 
The model can also be written in more compact form as 



(J)(B)z^ = a^ , 



where 



^ (|)(B) = (1 - (|)^B - - ... - , 

a polynomial in the operator B. Then, c() (B) represents a 
general non-seasonal AR operator. The reason these models 
are called autoregressive can be seen from recalling 
regression analysis; here the model relates the dependent 
variable to a set of explanatory variables 
...,z^_p), which are previous values of the time series. 
Therefore, the model is autoregressive. For example, when 
the value of the time series at time t depends only on the 
value at time t-1, the model becomes an autoregressive model 
of order 1, denoted AR(1) . Similarly, if the value in 
period t depends on the values in periods t-1, t-2 and t-3. 
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the model is AR(3) . The following equations are the mathe- 
matical representation of those two examples; 

AR(1) ; 

AE(3) : . 

2 . Moving Average Models 

In moving average (MA) models, the assumption is made 
that the current value of the time series at period t can 
be expressed as a linear combination of the previous fore- 

/N 

cast errors (or, residuals) , a^ = . The general 

equation for this model can be written 



^t ®l^t-l ®2^t-2 



0 a. 
q t-q 



/ 



which is called a moving average model of order q, where 
6g represents the scalar weighting coefficient for the q 
previous period. Like the AR model, the MA model can be 
written in more compact operator notation as follows : 



th 



z^ = 6(B)a^ , 



where 



9(B) = (1 - e^B - e^B^ - ... - e^B^) , 
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a polynomial in the operator B. The MA model implies that 
the analyst can gain valuable information for future 
predictions by consideration of the weighted sum of a 
number of previous forecast errors or residuals. The 
following are examples of MA models of order 1 and order 
3, MA(1) and MA(3) : 



MA(1) : = (1 - 9j_B)a^ = a^ - 

MAC3) : = (l-e^B-92B^-02B^) a^ 



^t " ®l^t-l " ®2^t-2 ~ ®3^t-3 



3 . Mixed and Integrated Autoregressive-Moving Average 
Models 

A natural extension of the autoregressive and moving 
average models is to construct a combination of the two. 

Such mixed processes are called autoregressive-moving average 
(ARMA) processes, of orders p and q, often written as 
ARMA(p,q) . Box and Jenkins [Ref. 4] notes that for many 
series encountered in practice, the inclusion of both AR 
and MA terms in a model results in fewer total parameters 
than would be required for either a pure AR or pure MA 
process. The ARMA(p,q) model may be written in the usual 
form as 









. + 4> z , + a. — 9, 3- ^ — 9^3 _ 

^p t-p t 1 t-1 2 t-2 



-9 a 
q t-q 
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in the more compact operator notation, the same model is 
written 




0 (B)a^ . 



As an example, an ARMA(2,2) model would be written as: 







or expanded, as: 



z 



t 




Recall that with differencing or other transformation 



techniques, a stationary series can be obtained from a 
homogeneous (seasonal or non-seasonal) non-stationary series. 
When the differencing technique is integrated directly into 
the ARMA model, the result is an integrated autoregressive- 
moving average (ARIMA) model. The model is often written 
ARIMA(p,d,q) , where p and q retain their previous meanings, 
and d represents the order of (non-seasonal) differencing. 

The model may also be generalized by including an additional 
constant term 9^ that will have the effect of adding a 
deterministic polynomial trend of order d. The value of 9^ 
may be found in general by the formula 
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/\ 



e 

o 



/V. 

= y G 



where 



P P 

G = [1 - I i^][l - I Sj] . 

i=l j=l 

The general ARIMA(p,d,q) model can be written compactly as 
(1) = 0^ + 6(B) a^ . 



An equivalent form for ARIMA(p,d,q) models that is 
sometimes seen is called undifferenced form: 



'f'(B) Zt " ®o ^t ' 



where 



'i'(B) = 4>(B) = (j)(B)(l - B)*^ . 

The "differenced" form, equation (1) , is the form usually 
seen. In the differenced form, the transfer function <}> (B) 
is assumed to be stable; that is, all the roots of the 
difference equation 4) (B) = 0 are outside the unit circle. 
Clearly, '{'(B) is not a transfer function of a stationary 
series when d > 0, since '{'(B) has d roots on the unit 
circle. 
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4 . Seasonal ARIMA Models 

The general class of ARIMA models previously dis- 
cussed can be further generalized to allow for the modeling 
of seasonal models with the addition of appropriate seasonal 
operators and parameters. The seasonal model is called 

an ARIMA model of order (p,d,q) x (P,D,Q) , where p,d, and 

s 

q retain their previous meanings, and P, D, Q and s are as 
defined in Section A of this chapter, page 17. The ARIMA 
(p,d,q) X (P,D,Q) model can be written as the product of 
nonseasonal and seasonal operators: 



As an example, consider the ARIMA (1,2,1) x (1,2,1) ^2 niodel: 



F. THE AUTOCORRELATION AND PARTIAL AUTOCORRELATION FUNCTIONS 
Mathematically, the autocorrelation at lag k, denoted 
Pj^, is defined as: 



The autocorrelation function of a time series describes the 
association (mutual dependence) among values of the same 



4)(B)$(B®)vMz^ 



9 +0 (B)0(B®)a^ . 

o ' ' ' ' t 




8c <l-8lB)(l-Sl,12®“>^t • 



- U) (Zt+k ■ 
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series taken at different time periods , the difference in 
time being referred to as the lag, denoted k. Autocorrela- 
tion of itself implies nothing about a change in one 
variable causing a change in another. However, the auto- 
correlations may provide important information about the 
structure of a data set and its pattern. In a set of com- 
pletely random data the theoretical autocorrelations among 
successive values will always be zero, whereas data values 
of strong seasonal or cyclical behavior will be highly 
autocorrelated. This pj^ can be estimated using the time- 
averaged sample autocorrelation, denoted r^^, defined as; 



N-k _ _ 

1 I (z - z) (z - z) 
N t=l ^ ^ ^ 



A plot of the autocorrelation function versus the lag k, 
called a correlogram, is very useful for the purpose of 
determining if a process is stationary and for identifying 
the appropriate model. 

Another function useful in identifying the appropriate- 
model for a given series is the partial autocorrelation 
function, often called the pauto function. Partial auto- 
correlations are analogous to autocorrelations in that they 
indicate the relationship of the values of a time series to 
various time-lagged values of the same series. However, 
they differ from autocorrelations in that they measure the 
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strength of the relationship between values of the series 
of various lags after the effects of other lags have been 
removed. In effect, they show the relative strength of the 
relationship that exists for varying time lags. For time 
series models of the types to be considered, the partial 
autocorrelation coefficients can be calculated several 
ways. One method of calculation, described by Box and 
Jenkins [Ref. 4, Sec. 3.2.5] begins with the following 
equation, satisfied by the autocorrelation function, where 
the are autoregressive parameters: 






* ♦k,k-l''j-k+l ♦kk'^j-k ' 3 = 



This leads to what are known as the Yule-Walker equations , 
which may be written 



1 


» — 1 

Q. 


p2 ... 


Pk-1 




^kl 




" Pi" 


Pi 

JL 


1 

• 

• 

• 


• • • 

* 

• 


Pk-2 

• 

• 

• 




'J>k2 


= 


P2 

• 

• 

• 


Pk-i 


• 

Pk-2 


• 

Pk-3 •• 


1 




<t> 

kk 




• 

Pk 



Solving these equations for k= 1,2,3, ..., successively, 
we obtain 



<P 



11 



r 
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P2 - Pi 
1 - P? 
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1 

P 



1 

1 



, etc . 



In general, for the determinant in the numerator has 

the same elements as that in the denominator, but with the 
last column replaced with The quantity a function 

of the lag k, is called the partial autocorrelation function. 
Other methods for calculating partial autocorrelations include 
successively fitting autoregressive processes of orders 
1, 2, ..., k by least squares and picking off estimates of 






11 ' '*'22 



, of the last coefficient fitted at each 



stage. Another method, discussed in Appendix A3. 2 of 
Ref. 4 and due to Durbin, generates estimates of an auto- 
regressive process of order k+1 recursively from estimates 
of autoregressive processed orders k and less. It is 
derived by observing recursive relationships from the 
Yule-Walker equations. The recursive formulas are: 
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‘*^k+l,j " ‘*^k+l,k+l'*’k,k-j+l ' 

j = 1,2, ..., k, and 



^k+l,k+l 

Another way of looking at partial autocorrelations 
is to consider a time series {z^} and the model 

^t " ^o ^l^t-1 ^2^t-2 \-l^t-(k-l) ' 

where- the b values are the least squares estimates of the 
linear regression coefficients (S's) in the model 



y ^ 

k+l j^'i^'^kj ^k+l-j 



^ - 



^t ^o ^l^t-1 ■*■•••+ ^k-l^t-(k-l) ®t • 



Let be the residual of after removing the linear 

~ /s 

effect of z. , , ..., z . . , , from z. , such that z. = z. - z. . 

t-1 t-k+1 t t t t 

Then, the partial autocorrelation of lag k, denoted 
is defined to be the simple autocorrelation of lag k for 
the adjusted series (z^,Z 2 , . 

The exact expression for partial autocorrelations of 
a moving average process is complicated, but an approximate 
duality between autocorrelations of an autoregressive 
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process and pautos of a moving average process implies 
that the pautos of a moving average process would have the 
same general behavior as the autocorrelations of an auto- 
regressive process of the same order. This approximate 
duality also implies that the pautos of an autoregressive 
process of order p will have the same general behavior as 
the autocorrelations of a moving average process of order 
p. Thus, one function can be examined and used to confirm 
the other, and may be of use in model identification, discussed 
in section G of this chapter. 

G. COMMENTS ON MODEL SELECTION EMPLOYING AUTOCORRELATIONS 

AND PARTIAL AUTOCORRELATIONS 

The purpose of this section is to provide a brief over- 
view of the methodology for selecting the proper model for 
a given time series , through both examination or prior 
knowledge of the data itself and examination of plots of 
the autocorrelation and partial autocorrelation functions. 

A not yet stationary series should first be differenced 
or transformed until it appears to be stationary; character- 
istic of the autocorrelation plot (correlogram) of a non- 
stationary series is a very slow damping out of the auto- 
correlations. When this property of the correlogram is 
observed, the analyst should difference the series until no 
further improvement is reached. For ARMA(p,q) models, the 
tentative identification of the model class, that is, the 
determination of p and q, can be done by examining the 
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sample autocorrelation and sample partial autocorrelation 
functions of the given time series with the theoretical auto 
and pauto functions of members of the general linear class. 
For most stationary time series, an adequate fit can be 
found in a model with p and q relatively small, say three 
or less. 

The following figure displays examples of autocorrela- 
tion and partial autocorrelation versus lag plots for several 
basic model types, described in Ref. 17. 

Figure 3. Examples of Auto and Pauto Plots 
Figure 3(a). Auto and Pauto vs. Lag for AR(1) Model 
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Figure 3(b). Auto and Pauto vs. Lag for AR(2) Model 
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Figure 3(c). Auto and Pauto vs. Lag for MA(1) Model 
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Figure 3(d). Auto and Pauto vs. Lag for MA(2) Model 
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Figure 3(e). Example Auto and Pauto vs. Lag for 
Mixed ARMA(1,1) Model 
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As can be seen from the figures, there is a wide variety 
of possibilities for appearance of the correlogram and 
partial autocorrelation plots for even the simplest of 
models. Sometimes the pattern of estimated autos and pautos 
can be quickly classified in terms of one of the simpler 
basic models; however, most real data will generate auto 
and pauto plots that resemble those of Figure 3 (e) , those 
of a mixed model. Here some free association may be re- 
quired to infer a pattern from the correlogram, or more 
than one pattern may be implied. However, most of the time 
a selection of something is possible. As is evident, this 
is a somewhat subjective process, where the quality of 
identification will improve with experience. 

The following table, adapted from Box and Jenkins 
[Ref. 4], summarizes the properties of autoregressive, 
moving average and mixed ARMA processes; an understanding 
of its contents will be helpful as thumbrules for initial 
model identification. 

It can be noted from examination of Table II that there 
exists what can be termed a "duality" relationship between 
the autocorrelations and partial autocorrelations of pure 
autoregressive and moving average processes. For example, 
the plot of autocorrelation versus lag for an AR(p) process 
would appear the same as the plot of partial autocorrela- 
tion versus lag for an MA(q) process, for p = q. Similarly, 
the plot of partial autocorrelation versus lag for the AR(p) 
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Autoregressive 

Process 


Moving Average 
Process 


Mixed 

Process 


Model in terms 
of previous 
2 ^ values 


(j) (B) 2^ = a^ 


6(B)a^ 


(|) (B) = 9 (B) 


Autocorrelation 

Function 


Infinite and tails 
off; conposed of 
dairped e:^nentials 
and/or daiiped sine 
waves 


finite; that 
is, there will 
be q non- 2 ero 
autocorrelations 


Infinite and 
tails off; 
cotposed of 
danped exponen- 
tials and/or 
danped sine waves 
after the first 
q-p lags 


Partial 

Autocorrelation 

Function 


finite; that is, 
there will be p 
non- 2 ero partial 
autocorrelations 


Infinite and 
tails off; 
daninated by 
danped exponen- 
tials and/or 
sine waves 


Infinite and 
tails off; 
conposed of 
danped exponen- 
tials and/or 
sine waves after 
the first p-q 
lags 



Table II. Properties of Basic ARMA Processes 



process would appear the same as the plot of autocorrelation 
versus lag for an MA(q) process, where p = q. 

For seasonal series, where differencing has been unable 
to remove all the nonstationary seasonal behavior, the auto- 
correlation and partial autocorrelation plots will generally 
exhibit a "spike" at lags equal to integer multiples of the 
period of seasonality. For example, if the data are seasonal 
with a period of 12, then the autocorrelations 
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r, r_-, ... and the partial autocorrelations 
Iz Zh Jo 

d),_ d)«. d)-- ... would be amplified. This would 

XZ/J.Z Z ^ f Z ^ JO/Ju 

indicate a need for one or more seasonal parameters to be 
included in the model. The following figure is an illus- 
tration of such an auto/pauto plot. 

Figure 4. Example Auto and Pauto versus lag for mixed 
ARMA(p,q), with seasonality remaining of 
period s = 6. 
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H. ESTIMATION OF MODEL PARAMETERS 

Suppose a time series has been tentatively identified 
as an ARIMA(p,d,q) model: 

(}) (B) = e + e (B) a^ . 

t o t 

There are then p+q+2 unknown parameters of the model to be 

2 

estimated, ((}>.,(})_, (|) ,9 ,9^,9-, ..., 9^, a ). The 

1, Z pO-LJ < 53 . 
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Box-Jenkins procedure separates this estimation into two 
parts. First, estimates are obtained for the AR and MA 

^ /N 

P3.3T3.in0 1 03T S ^ f ^ f m m m f 3nd 9 “• •••/ 

2 

then estimates are calculated for 9 and a , which are func- 

O 3 

/N /V 

tions of (p and 9 . 

The usual procedure is to select those parameter values 

4> and 9 that minimize the sum of squared model errors 

(residuals). Let = 9^/(1 “ “ ‘^^2 ” *’* ~ *^p^ 

= V^(z^ - y ). Then, it can be shown that = E[W, ] 
t t w w t: 

and that the model can be rewritten as 

4>(B) (W^ - y^) = 9(B) a^ , 



or 



a^ = 9‘^(B) ({)(B) (W^ - y^) . 



Now, let W, the sample mean, be the estimate of y^ (if 
d > 0, then y^ is usually zero) . Also, let 



a^ = 9"^ (B) $ (B) (W^ - W) . 



Set 

N 

S(±rh = I , 

t=l 

where N = n-d is the length of the differenced series. The 
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objective is to select those parameters ^ and £ such that 

/N. 

S(^,^) is minimized over all values of the (p+q) -dimensioned 

parameter space. Since the equation in S is nonlinear 

in the parameters, iterative search methods are generally 

used to determine optimal model parameters. Having found 

2 

these optimal parameters and the minimum S, o can be 

3 . 

estimated using 



Sf = S($,0)/N-p-q . 

Chapter III contains a more detailed description of a 
nonlinear algorithm to estimate parameters of a generalized 
seasonal or nonseasonal ARIMA model. 

I. DIAGNOSTIC CHECKING OF MODELS 

After the model has been tentatively identified and 
parameter estimates have been calculated, the next task is 
to test whether or not the original model specification 
was correct and that the model itself is adequate in fore- 
casting power. The process of testing the model can take 
many forms, but will usually include at least the following 
two steps: 

1. Generate a simulated series from the estimated 
model and compare the simulated series and its 
auto and pauto functions with the original series 
and its respective auto and pauto functions. This 
comparison is essentially subjective. 
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2. Calculate the residuals of the estimated model, the 

/s 

a^'s, and compare the properties of the residuals 
with those assumed for the shock terms of the actual- 
process. The residuals should be normally distri- 
buted and uncorrelated with each other; that is, there 
should be no discernable structure in the residuals. 
There are many quantitative statistical tests and 
data analytic tools that can be applied to the 
residuals to test hypotheses of normality and zero 
autocorrelation . 

A plot of the autocorrelation and partial autocorrelation 
functions of the residuals will provide not only a test of 
whether or not the residuals are uncorrelated, but, if they 
are correlated, the plots can be used to suggest improve- 
ments to the model. For example, suppose the model was 
tentatively specified as the ARMA(1,1) model below: 

(1 - 0.5B) (z^ - 2) = (1 + 0.7B) a^ . 

Suppose also that the autos and pautos of the model residuals 

suggested the model (1 - 0.3B) a^ = u^, where the 

2 

are whxte noise, uncorrelated with variance a . Then, 

u 

these two models can be combined into the ARMA(2,1) model: 

(1 - 0.3B) (1 - 0.5B) (z^ - 2) = (1 + 0.7B) u^ . 

This should be a refinement over the original model. 
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J. FORECASTING 



The objective in forecasting is to predict future 
values of the time series with as little error as possible. 
The criterion most often used for forecasting is to calcu- 
late that forecast which minimizes the expected mean square 

/N 

forecast error. Therefore, if z^(iJ.) denotes the forecast 
for lead time I from time origin t of the value 

/N 

the objective is to find z^{l) such that the objective 
function 



Et (z 



t+£ 



2 



is a minimum. This forecast is given by taking z^(i) as 
the conditional expectation of given z^, ..., z^: 



z^U) 



= E [z 



t+il 



't-1' 



^ 1 ^ 



The forecast can be easily generated recursively from the 
mathematical model, utilizing the facts that 



“ ^t-j ' 



for j = 0,1,2... (where t is the current time), and 



E[z^_i_j I z^, . . . ,z^] = z^(j), for j = 1,2, 
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and 



0 for times > t 

E[a^] = 

a^ for times £ t 

For excimple/ consider the Box- Jenkins model 
(1 - 0.5B + 0.6B^) = (1 + 0.3B) a^ 



which could be written in the expanded form as 



z 



t 



0.5z 



t-1 



0-6^t-2 * \ ®t-l 



Now, assume it is known that Zq^qq is 1.4, z^^ is 1.0, and 
the calculation has been made to obtain the residual 
a^QQ = 0.2. Then the forecasts of z^q^, ^-^ q 2 ^103 

made from an origin of t = 100 can be calculated from the 
model as follows: 

first, for £ = 1, set 



^100 



( 1 ) 



E [z 



101 ' 100' 99 




/ 



^100 



( 1 ) 



E to. 5z 



100 



- 0 . 6z 



99 



+ a 



101 






and 
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(0.5) (1.4) 



(0.6) (1.0) + 0 + (0.3) (0.2) 



^100 



then our forecast 



^qq(I) = . 0.7 - 0.6 + 0.06 = 0.16 ; 



Next, for 1=2, set 



^100 



^^^102 I ^ 101 '^ 100 '^ 99 ' * * * '^ 1 ^ ' 



^ 100 ^^^ ^^°‘^^101 ” °*^^100 ^102 ^*^^101 



^100'^99' * * * '^1^ ' 



z^qq(2) = (0.5) (z^qq(D) - (0.6)(Zj_qq) + 0 + 0 , 



Zioo(2) = (0.5) (0.16) - (0.6) (1.4) = -0.76 ; 



similarly, for i = 3, set 



^100^^^ ^^^103l^l02'^101'^100,^99' * * "^1^ ' 



^ 100 ^^^ E [ 0 . 5 Zj _ q 2 ■ ®*®^101 ^103 °*^^102 



^100'^99' • • • '^1^ 



Zioo(3) = (0.5) (z^QQ (2) ) - (0.6 ) (z^qq(1) ) + 0 + 0 , 
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and then 



ZioqO) = (0.5) (-0.76) - (0.6) (0.16) = -0.48 . 

Clearly, this process can be continued into the future as 
long as desired, recognizing that the expected forecast 
accuracy will decrease as I increases. 

Let 



et(^) = 



denote the forecast error for lead time I beyond the forecast 
time origin t. It can be shown that e^(i) is given by 
the relation 



e^(^) = 



^t+il 



^l^t+Ji-1 



'^5,-l^t+l 



where the scalar weights are determined from the equation 



t|;(B) = $ ^(B) (1 - B) ^ 0(B) 



using the known parameter values. The variance of the 
forecast error is given by 

Si-1 

E[(e^(il)] = (1 + ipj + ... + (t^) = {1 + I . 

j=l 
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From these relationships, a confidence interval of n 
standard deviations about a forecast of lead time I would 
be given by: 



Cn(Zt ) 



z^a) 






n [1 + 



1-1 



I 

j=l 



1/2 



a 

a 



The value of n could also be though of as the percentage 
point for the desired level of confidence using the Normal 
distribution [Ref. 4]. It can be noted from the definition 
of above that the one-step-ahead forecast error, 

e^(l) is simply the value that is. 



=^t+i - 



®t'^> = ^t+i 



This is an explanation for the use of the term "residual" 
to refer to the white noise or shock terms. Additionally, 
from the foregoing it is evident that the forecast error 
variance is a non-decreasing function of the length of 
the forecast lead time therefore, the confidence bands 
must become wider as the forecast lead time increases. 
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III. GENERALIZED ARII4A MODEL PARAMETER ESTIMATION 



The purpose of this chapter is to describe the compu- 
tational details of the algorithm WMARQRDT, employed with 
appropriate subroutines to estimate Box-Jenkins parameters 
of a generalized seasonal (or non-seasonal) ARIMA model. 
This generalized model can be written parsimoniously as 
follows ; 



(j) (B)<&(B®)W^ = 9^+9 (B)0(B®)a^ , 



where 



W 



t 



is V*^V°z 



d t 



r 



and 



z^ is the transformed time series (for example, 

z' = In z. or z' = . 

t t t t 

The model employs a Marquardt-type (see Ref. 13] non- 
linear least squares algorithm for determination of the 
model parameters, searching in the parameter space for the 
set of parameters which minimizes the sum of the squared model 
residuals . 
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A. MATHEMATICAL DISCUSSION OF THE MARQUARDT ALGORITHM 

Reference 13 describes Marquardt's methodology for 
development of a general algorithm for least-squares esti- 
mation of nonlinear parameters. Most algorithms for least- 
squares estimation of nonlinear parameters have been cen- 
tered about either a pure linear iterative model based on 
a Taylor series expansion or using some form of the method 
of steepest descent (maximum negative gradient) . Since 
both methods have severe potential pitfalls (Taylor series 
due to divergence of the iterates due to an unfortunate 
choice of initial values, the steepest descent due to slow 
convergence after the first few iterations) , a method termed 
"Maximum neighborhood" method was developed by Marquardt 
which, in effect, performs an optimvim interpolation between 
the Taylor series method and the method of steepest descent , 
the interpolation being based upon the maximum neighborhood 
in which the truncated Taylor series gives an adequate 
representation of the nonlinear model. 

As discussed in Ref. 14, the idea of Marquardt's method 
can be briefly explained as follows. Assume the algorithm 
begins with a vector of initial parameters *^o ' ®o ' ®o^ * 

If the method of steepest descent is applied, a certain 
vector direction h^, where the subscript g represents the 
gradient (which gives the direction in which the rate of 
change is the greatest) , is obtained for movement away from 
the initial point values in the parameter space. Due to 
potential nonlinearities in the sum-of-squared-residuals 
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surface, ^ *^o' parameter space, may 

be the best local direction in which to move toward opti- 
mality (minimal sum of squared model residuals) , but may 
not be the best direction. However, the best direction 
must be within 90 degrees [Ref. 8] of h , or else the 
value of S (({),$, 9,0) will increase locally. Now, the 
linearization, or Taylor series method, .may result in a 
different correction vector, h^. Marquardt found experi- 
mentally that for a number of practical problems, the angle, 
say C/ between the vectors h^ and h^, fell in the range 
80 degrees < ^ < 90 degrees. In other words, the two 
directions were nearly normal in most cases. The Marquardt 
algorithm then, provides a method for interpolation between 
the vectors h^ and h^, employing a sutiable step size, which, 
in general, converges more quickly than either of the two 
"parent" methods alone. The remainder of this section des- 
cribes the mechanics of the algorithm, as applied to the 
estimation of Box-Jenkins parameters. 

B. MODEL INPUT INFORMATION 

The following is the required input information to 
the WMARQRDT parameter estimation model: 

{W^} the transformed time series (for example, 

= In z^, or = /z^ 

N the length of the time series 

s the length of the seasonal period (for no 

seasonality, s = 1) 
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p 

p 

Q 

q 

NDIFNS 



NDIFS 







the number of seasonal AR parameters, P ^ 0 
the number of non-seasonal AR parameters, 

p ^ 0 

the number of seasonal MA parameters, Q ^ 0 

the number of non-seasonal MA parameters, 
q ^ 0 

the number of non-seasonal differences to 
be taken 

the number of seasonal differences to be taken 



$p initial estimates of seasonal AR 

' parameters 






initial estimates of non-seasonal 
AR parameters 



0 



^®l,s'®2,s'*-’ '®Q,s^ 



initial estimates of seasonal 
MA parameters 



/N. 

6 = ( 9 ^, 02 ^ 




initial estimate of non-seasonal 
MA parameters 



C. ALGORITHM DESCRIPTION 

Initially, the model accepts inputs of the time series 
length and the (transformed, but undifferenced) time series 
itself from the user's disc, where it is stored in File 
FT02F001. The remainder of the input data is then entered 
interactively on the terminal as prompted by the model . The 
model then forms the initial parameter estimates into vectors 
of lengths p, P, q, and Q, and employs them to calculate 
an initial model residual sum of squares. 
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Prior to initiating the sum of squares calculation pro- 
cess, the model unravels the polynomial of AR and MA operators. 



and forms a vector of coefficients for both the AR and MA 
sides of the model. For example, consider the relatively 
simple (2,1) x (2,1) ^2 -ARIMA model, written 



This is represented in the computer as four vectors : two 

vectors, tt and y, composed of the expansion AR and MA 
coefficients respectively, and two vectors composed of the 
indices of the time series terms associated with the coeffi- 
cients. For the model in the example above, these vectors 
would be written as follows : 





Unraveling, the model becomes first: 



(l-<j)3_B-(|)2B^) (W^-$ 



w 

l,12^t-12 



(1-9^B-02B^) (a^-0 



l,12^t-12 



and then takes the form: 





TT 




INDEX 



AR 



(0, 12, 1, 13, 2, 14) 



Y 
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INDEX. 



( 0 , 12 , 1 , 13 , 2 , 14 ) . 



MA 

This format for unraveling of a generalized product of 
polynomials of seasonal and non-seasonal parameters permits 
simplified calculations in the determination of the sum of 
squared residuals, since the arithmetic need only be per- 
formed for the non-zero parameter values, those contained 
in the n and y vectors. In the actual calculations, per- 
formed in FORTRAN, the index vector origin is translated 
to the right by the amount p + Ps + q + Qs, called lADDIT 
in the model, in order to prevent obtaining zero or negative 
subscripts . 

The calculation of the residual sum of squares is a 
three stage process. First, the model calculates ten 
values of e^, the "forward” white noise terms, using the 
following relationship: 

p+Ps q+Qs 

e. = W. - y ir.W, . + y Y-e. .. 

t t ^ 1 t-i ^ 3 t+3 

i=2 j=2 

The assumption is made that equal 

to zero. Having found the desired values for (e^Q ,e^ , . . . , 0 ^^) , 
we can now backcast the necessary values of W^, in order 
to enable us to calculate the required white noise terms 
a^, as follows: 
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w. 



P+Ps 

i=2 



g+Qs 

I Y-;® 
j = 2 



j t+j 



solving for " '^1- (p+sP+q+sQ)^ ' 

W^'s backcast, we can proceed to calculate the estimated 
white noise terms (residuals) using this relationship; 



= W - 



p+Ps 
I 7T^ W 
i=2 



i t-i 



q+Qs 

+ I Y 
j=2 



j^t-j 



The values of a^ are calculated for t = 1, . . . , N. Once the 
a^'s are calculated, it is a simple matter to calculate the 
sum of squares, using 



N 

S(4),$,9,0) = I • 

t=l 



Having calculated the residual sum of squares for the 
initial parameter estimates, we form the parameters into 

/V 

a vector 3, such that 



6 ~ B 2' ' • • f f 



where k = p+P+q+Q ; that is. 



3 



(<!> 



1 ' 



''^p''^l,s' 






P,S 



/ 



e 



1' 



/9 f 

q d,s' 



,e 



Q,s 



) . 
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The subroutine PARSH then calculates the derivatives 



X t 

i,t 36 ' 

i 

over all values of t = 1/.../N and i = l,...,k. Using the 
residuals calculated earlier, the derivatives are 

estimated numerically using a perturbation of each parameter 
of amount 6 , such that 



- {a^(6^, . . . ,Bj^)-a^(6j_/ . . . /B^+6 






Then with {a } and x. supplied for the current parameter 

Ti X / u 

values, the following quantities are formed: 

1. The kxk matrix A = j ^ ' where 



N 



A. . 
ID 



t=l 



2. the vector G = {g^^} , where 



N 



9i = I 
t=l 



3. the scaling quantities = /A^^ 



58 



Then, the modified (scaled and constrained) linearized 
equations 



. [A*] [h*] = [g*] 

are formed, where 

A. .* = A. ./D.D. , i 7^ j , 

A..* = A..*+ 5, 5=.01, 

11 11 ^ 



and 



* — 



g./D. 
^1 1 



It can be noted here that the [A^^*] matrix is in fact a 
correlation matrix for the model parameters. Those equa- 
tions are then solved for h*, which is scaled back to give 
the parameter correction vector h^ , where 



hj = h.VDj 



Then, the new parameter values are formed from the h vector, 
such that 



»NEW = ®OLD ' 
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4 






i 





and the sum of squared residuals for the set of 

NEW 

parameters is calculated. 

NOW/ if S(3„_J < the parameter corrections h 

NhW ULiJJ 

are tested. If all elements of the h vector are smaller 
than some epsilon, say .0001, then convergence is assiimed 
and the k x k matrix [A*] ^ is used to calculate the covari- 
ance matrix of the estimates; otherwise, is reset to 

^NEW' constant ? is reduced by a factor F, say .1, and 

computation returns to the calculation of a new set of 
derivatives . 

However, if 3(3^^^^) > S(3 qj^j^), the constraint parameter 
5 is increased by the factor F, and computation is 
returned to the stage where the A matrix is formed. In 
all but unusual cases, a minimal sum of squares will be 
found. However, an upper bound can be placed on and 
when it is exceeded, the search is terminated. When con- 
vergence is reached, the desired output information is 
calculated. 

D. MODEL OUTPUT AND DIAGNOSTIC CHECKS 
The output includes; 

1. 3 = (32/* the least-squares estimates of 



' * * • ' 




9 



1 



/ • • • f 



0 



q 



0 



1 



r • • • r 



0 



Q 
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2. a‘ 



the residual variance. 



N-p-q-P-Q 



S((J),$,0,0) 



3. V 



the covariance matrix of the parameter 
estimates, formed from 

nr ^ 

V = {V. . } = (A A) a 
13 a 



4. s. 



the standard errors of the parameter 
estimates. 



Si = , i = 1, . . . ,p+q+P+Q 



5. 



the correlation matrix, obtained from 



R.. = V../ /V..V.. 
13 13 11 33 



6. e 



the overall constant term, where 

0 = y G , for 

o 



P P 

G = (1 - I (1 - I ^j) 

i=l j=l 



and 



N 



I 



t=l 



Finally, diagnostic checks on the model are performed. 
This includes calculation of residual autocorrelations 
(k) , obtained from these formulas: 



r^^(k) = c^^ (k) /c-'^ (0) , 
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where 



N-k 

“ F I (it - <®t+k • ' 

t=l 



N 




t=l 



and k goes from 1 to the maximum desired lag K, usually 
40. Additionally/ a chi-square statistic is calculated 
from 



X 



2 

(V) 



K 

= N I , 

k=l 



and is compared with a chi-square distribution with 
v=K-p-q-P-Q degrees of freedom. 

Examples of input methodology and appearance of output 
for this model are contained in later chapters of this 
paper . 



62 



IV. DESCRIPTION OF THE TIME SERIES EDITOR 



In this chapter, descriptions are given for each of 
those programs in the Time Series Editor which interact 
with the user. There are fifteen separate program modules 
currently included in the Editor; there are other programs 
utilized by the Editor that are completely transparent 
to the user, and are not described here. However, listings 
for all programs contained in the Editor are provided in 
Appendix D to this report. 

No attempt is made in this chapter to describe the 
actual mathematical calculations or algorithms that are 
performed by the programs. Rather, the objective is to 
provide general descriptions of what each program can do 
for the user, and how the user interacts with the programs. 
The programs described in this chapter are, in the order 
presented; TIMESER EXEC, ZFORMAT, CMSWORK, TRANS, DIFF, 

PLOT, AUTO, ESTIMATE, YESTSEAS, WMARQRDT, XSUMSQ, FORECAST, 
ROOTS, SIMULATE and GENERATE, and HELP. 

A. THE EXECUTIVE PROGRAM 

The heart of the Time Series Editor is a master program 
called TIMESER EXEC which provides file control for all of 
the other program modules, controls input and output, super- 
vises the necessary CP/CMS protocol and provides instructions 
interactively to the user concerning the contents of the 
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Time Series Editor and how each program can be used. 

The TIMESER EXEC program is written in the IBM-360 CP/CMS 
Executive language. It is the only program in the package 
not written in FORTRAN; consequently, it should be the only 
program that would need modification if the Editor were to 
be adapted to another FORTRAN-capable time-sharing system. 

After the user has logged into CP/CMS and linked to 
the disc space containing the Time Series Editor (instruc- 
tions for this procedure are provided in detail in the 
User's Guide to the Time Series Editor, included as Appendix 
A to this report) , the entire Editor package is made avail- 
able to the user by the command TIMESER. On entry of this 
comand, the executive routine, TIMESER EXEC, can provide 
a guided tour through the Editor. It will briefly des- 
cribe what the Editor can provide and asks the user what 
option he wants to use. On the basis of the user's response, 
the Editor then advises the user what input data is required 
and how it should be entered, either through a data file 
entered offline or parametric data entered interactively 
at the terminal. When the user selects an option for exe- 
cution, the TIMESER EXEC routine loads the appropriate pro- 
gram package and automatically manipulates the required 
input and output data files. For the more experienced user 
who does not require detailed user instructions, there is 
a shortened version of the TIMESER EXEC. Upon entry into 
the shortened version, the user is immediately asked which 
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option he desires. This version of the Editor is entered 
by adding the argument "s" when logging into the Editor; 
that is, type TIMESER S when logging in. The TIMESER S 
(shortened) version of the Editor provides exactly the 
same program options to the user as the longer TIMESER 
version. 

B. DATA AND FILE MANAGEMENT 

Whenever data is required, the user is prompted at 
the terminal by either the TIMESER EXEC or the program 
module being executed. In most cases, the necessary user 
response is a short alphanumeric character input during 
execution using the terminal keyboard. However, when the 
time series data itself is entered, the user may enter the 
data by card deck offline. Similarly, most output is pro- 
vided right at the user's terminal. Output such as listings 
of transformed series or plots are sent to the offline 
printer to conserve time and provide hard copies of the 
results. Detailed descriptions of the input and output 
requirements of each program module are provided in this 
chapter, and also in abbreviated tabular form in Appendix 
A, the User's Guide. However, there are some general prin- 
ciples of data management that apply to all programs . 

These are described in this section. 

1. Offline Data Entry 

The Time Series Editor requires that the user's 
time reside in FILE FT02F001. The data can be read offline 
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via cards and directed to FILE FT02F001. An example of the 
required format for the input card deck for the Time Series 
Editor is provided in Figure 5. It is noted here that the 
data format expected by all programs is FORMAT (5F15 . 6) , 
preceded by a card containing the length of the time series 
in FORMAT (13) . 

2 . ZFORMAT Program 

In many cases, the user will want to analyze time 
series data that are already available ina format other 
than that required by the Editor. To spare the user the 
tedious task of retyping his deck in the required format, 
the ZFORMAT. program was written. This program converts a 
time series data file (FILE FT03001) written in any FORTRAN 
format into the proper FORMAT (5F15 . 6) onto FILE FT02F001. 
The length of the series is entered as an input on the 
terminal. The data are then ready to use in the Editor, 
with the original file left in FILE FT03F001. Figure 6 
shows a sample input deck for such non-standard data. 

Figure 7 gives a sample session where ZFORMAT is used. 
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Figure 5. Card deck arrangement for OFFLINE READ of time series input 
data in FORMAT (5F16. 5) 



I 



data in present format; 
here, it is 2F10.2 
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Figure 6. Card deck arrangement for OFFLINE READ of time series input 
data in any FORMAT 



timeser s 



ENTER LETTER FOR OPTION YOU WANT, 
z 

IS YOUR DATA IN FILE FT03F001? 

y 

EXECUTION BEGINS... 

IS YOUR TIME SERIES DATA NOW IN FILE FT03F001? 

y 



ENTER THE LENGTH OF YOUR TIME SERIES VIA FORMAT 13. 

010 

NOW ENTER THE FORTRAN FORMAT FOR YOUR TIME SERIES DATA, 
INCLUDING PARENTHESES; FOR EXAMPLE, TYPE: (5F15.7). 

flO. 2 

YOUR FORMAT IS: 

FORMATFIO. 2 
IS THIS CORRECT? 



NOW ENTER THE FORTRAN FORMAT FOR YOUR TIME SERIES DATA, 
INCLUDING PARENTHESES; FOR EXAMPLE, TYPE: (5F15.7). 

(flO.2) 

YOUR FORMAT IS: 

FORMAT (FlO. 2) 

IS THIS CORRECT? 

y 



YOUR DATA IS NOW TRANSFORMED INTO THE PROPER FORMAT 
FOR USE IN THE TIME SERIES EDITOR , LOCATED IN FILE 
FT02F001. YOUR ORIGINAL DATA IS IN FILE FT03F001. 

DO YOU WANT TO GO AGAIN? 
n 

CONTROL RETURNED TO CMS 
R; 



Figure 7 . 



Sample user session with ZFORMAT program 



3. Interactive Keyboard Input 



a. Numeric Input 

When the user is prompted at the terminal to 
enter such numerical values as the number of AR parameters, 
estimates of parameters or the number of differences to be 
taken, he must enter these values before program execution 
can continue. The user should enter the data according to 
standard FORTRAN practice. For example, integer data should 
be entered (without a decimal point) for counts, such as 
the length of a seasonal period, and for names beginning 
with the letters I through N; floating-point data, such as 
initial estimates of parameters, should be entered with a 
decimal point. Because a typed decimal point overrides a 
floating-point format, it is not necessary for the user 
to concern himself with the format for floating-point input 
data. However, care must be taken when entering integer 
data because it must be right- justified in its format field. 
The user is advised in each case when the format is other 
than II. For example, suppose the user desires to perform 
analysis on a series having only 90 observations. If the 
program that he is executing requires the length of the time 
series, the user will receive the following request; 



ENTER LENGTH OF THE TIME SERIES VIA FORMAT 13 



The user would then enter: 



column 123 

b9Q (where b represents a blank space) . 

If the blank were omitted, the program would read the 
series length as 900, and problems in execution would 
occur . 

b. Alphabetic Input 

The Editor often asks the user to respond with 
alphabetic input. Generally, this is in the form of the 
response to a question requiring a yes or no answer, or a 
response to the question of which program to execute. The 
editor has been programmed to read only the first letter 
of such responses. Therefore, the user need only enter the 
first letter of the response for each such inquiry. For 
example, the user should enter Y for yes, N for no, P for 
the PLOT option, X for the XSUMSQ option, W for the 
WMARQRDT option, etc. The other alphabetic input is 
generally entry of titles to plots; these can be any 
combination of numeric and alphabetic characters, as long 
as they do not exceed 72 columns in length. For example, a 
suitable title would be: 

PLOT OF AUTOS AND PAUTOS FOR SERIES C DATA (2 ORDINARY DIFFERENCES) . 

4 . Output 

Most of the results are written out at the user's 
terminal. However, in some cases the output is printed 
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offline to conserve time. Such results as plots and 
transformed time series are written onto various files 
(FT03F001, FT09F001, FT08F001) , and the TIMESER EXEC pro- 
gram causes them to be automatically printed offline under 
the user's identification number (USERID). The plots and 
files can also be printed out at the terminal using the 
usual CMS commands in the CMS mode, outside the Editor, or 
in the CMSWORK program (see next section) . However, the 
user should be aware that plotting graphs at the terminal 
is generally quite a slow process. 

5 . CMSWORK Program 

This program permits the user to enter the CMS 
environment to perform routine CMS functions without leaving 
the Time Series Editor. Commonly used CMS commands include 
file name alteration, file erasure, file offline printing 
or punching, obtaining a disc status, and listing all disc- 
resident files. However, essentially all CMS commands can 
be issued while in CMSWORK. Figure 8 provides an example 
of the use of the CMSWORK program. 

C . TRANS PROGRAM 

The TRANS program takes a given time series in File 
F102F001 and performs a transformation upon each data ele- 
ment. The options include a change of scale, a power trans- 
formation, a shift of the origin, and a natural logarithm 
transformation. The original data is written onto a file 
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timeser s 



ENTER LETTER FOR OPTION YOU WANT, 
c 

ENTER DESIRED CP/CiMS COMMANDS, ONE PER LINE. 

WHEN FINISHED TYPE; &GOTO -QUES 
stat 

P (191); 30 FILES; 263 REC IN USE, 33 LEFT (of 296), 39% FULL (2 CYL) 
listf * ft02f002 
FILE NOT FOUND 
17.14.13 LISTF * FT02F002 
! ! ! E(00002) ! ! ! 
listf * ft02f001 

FILENAME FILETYPE MODE NO. REC. DATE 



SERG 


FT02F001 


PI 


3 


9/16 


SERC 


FT02F001 


PI 


5 


9/16 


LNSERG 


FT02F001 


PI 


3 


9/16 


FILE 


FT02F001 


PI 


3 


9/16 



cp q f 

FILES:- NO RDR, NO PRT , NO PUN 
erase file ftOSfOOl 
offline print file ft02f001 
offline print serg ft02f001 
offline print sere ft02f001 
erase file ft02f001 
stat 

P (191); 28 FILES; 240 REC IN USE, 56 LEFT (of 296), 81% FULL (2 CYL) 
&goto -ques 

DO YOU WANT TO GO AGAIN? 
n 

CONTROL RETURNED TO CMS 

R; 



Figure 8 



Sample user session with CMSWORK program 



called DATA FT02F001, while the transformed series is 
written onto FILE FT02F001. The parameters of the trans- 
form itself are written onto FILE FT07F001. The TRANS 
program is self-contained, using no other programs. 

D. DIFF PROGRAM 

The DIFF program takes a given time series in FILE 
FT02F001, and performs non-seasonal and/or seasonal differ- 
encing operations on the series in order to help achieve 
stationarity for proper modeling. The user must indicate 
if seasonal differences are required; if so, the length of 
the season is input. The user is then required to enter 
the number of seasonal and nonseasonal differences. After 
execution of DIFF, the original data are in FILE FT02F001, 
and the differenced data are written onto FILE FT03F001. 

To use the differenced data in any other program, CMSWORK 
program must be entered and the files altered to save the 
original data (if desired) and to transfer the differenced 
data to FILE FT02F001. The CMS commands to perform this 
shift of files would look like: 

ALTER FILE FT02F001 PI SAVEDATA FT02F001 PI 
R; 

ALTER FILE FT03F001 PI FILE FT02F001 PI 
R; 

The DIFF program uses the IMSL subroutine called FTRDIF. 
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E. PLOT PROGRAM 



The PLOT program plots any given time series which 
resides in FILE FT02F001. Other than the time series, 
which must be on the user's disc or entered offline, the 
program requires only that an identification title for the 
plot itself be entered during execution. The plot is auto- 
matically printed offline. The PLOT program uses a modified 
PLOTP subroutine called PLOTS from the IBM Scientific 
Subroutine Package Library (SSPLIB) . 

F . AUTO PROGRAM 

The AUTO program takes a given time series in FILE 
FT02F001 and calculates summary statistics of utility in 
Box-Jenkins modeling. The statistics include the sample 
mean, variance, autocorrelations and partial autocorrela- 
tions for lags one though 40. Part of those statistics 
are printed out at the user's terminal, while plots of the 
autocorrelations and partial autocorrelations versus lag 
are printed offline. These plots will also be written onto 
the user's disc in FILE FT08F001, and may be plotted at 
the terminal in CMS outside the Editor or in the CMSWORK 
program; recall that this is a time-consuming process. The 
AUTO program provides the user with essential information 
about stationarity , seasonality and model identification. 
With the Box-Jenkins procedures, the second moments (auto- 
correlations and partial autocorrelations) are valuable 
tools in tentative model identification. Additionally, AUTO 
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is useful in the diagnostic checkout phase of model 
building, when examining the model residuals for structure. 

The AUTO program uses the IMSL subroutine FTAUTO. 

G. ESTIMATE PROGRAM 

After the user has tentatively identified a model for 
a series through analysis of the series plot, autocorrela- 
tions and partial autocorrelations, ESTIMATE should be 
executed to calculate maximum likelihood estimates of the 
model parameters. ESTIMATE is written to calculate parameters 
for non-seasonal models only; for seasonal models, the 
programs YESTSEAS and WMARQRDT should be used for parameter 
estimation. The ESTIMATE program requires that the series 
to be modeled (already transformed and differenced to 
achieve stationarity) reside in FILE FT02F001; additionally, 
the program will instruct the user to input the number of 
MA and AR parameters in the model, as well as the number of 
differences taken. The general model for which ESTIMATE 
calculates parameters is; 

+ 0(B) a^ 

where W' = V z' and z! is the transformed value of z. . 

t t t t 

ESTIMATE calculates estimates of the AR parameters for the 
undifferenced form of the model, the estimated MA parameters 
and MA constant term, the residual variance, autocorrelations 
and partial autocorrelations of the residuals, plots of 
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these autocorrelations and partial autocorrelations, and 
a Chi-square goodness-of-f it statistic for the model. The 
original series to be modeled remains in FILE FT02F001 
and the residuals of the model are stored in FILE FT02F001. 
ESTIMATE uses the IMSL subroutine FTMAXL (a nonlinear 
gradient search algorithm) to estimate the model parameters. 

H. YESTSEAS PROGRAM 

The YESTSEAS program is used to calculate initial esti- 
mates of autoregressive and moving average parameters for 
a generalized seasonal Box- Jenkins model. These initial 
estimates are used as inputs to the program WMARQRDT, which 
then refines them to deteimiine the parameter estimates that 
minimize the sum of the squared residuals. The general 
model assumed in YESTSEAS is 

(}){B) c&(B®)V^VgZ^ = 0 {B) 0 {B^)a^ . 

The YESTSEAS program requires that the seasonal time series 
for which initial parameter estimates are required is located 
in FILE FT02F001; it should be transformed, if desired, but 
generally not differenced, since YESTSEAS performs the 
required (seasonal and/or non-seasonal) differencing. 

Terminal inputs required by the program include the number 
of sdasonal and non-seasonal differences to be taken (if 
data is already differenced, enter zero for these values), 
the length of the seasonal period and the number of each 
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t 

type of parameter to be estimated (that is, the number of 
each of the seasonal and non-seasonal AR and MA parameters) . 
The output of the model, initial estimates for the requested 
parameters, is printed out at the terminal. YESTSEAS uses 
the IMSL subroutines FTRDIF, FTMAXL, FTAUTO and FTARPS, 
as well as a modified version of FTMAPS. 

I. WMARQRDT PROGRAM 

The WMARQRDT program uses the initial parameter esti- 
mates calculated by YESTSEAS or other means as starting points 
for calculation of parameters of a generalized seasonal 
Box- Jenkins model which minimize the sum of the squared 
residuals. The general seasonal model assumed by WMARQRDT 
can be written: 



4) (B)f (B^)V^V°z = 0 +0(B)0(B^)a. 

s t o t 

Details of the Marquardt-type non-linear least squares 
algorithm employed in this program are provided in Chapter 
III of this report. WMARQRDT requires that the series to be 
modeled reside in FILE FT02F001; it should be transformed, 
if desired, but generally not yet differenced, since the 
program performs differencing. Inputs at the terminal 
include number of seasonal and non-seasonal differences, 
length of the seasonal period, number and type of parameters 
desired, and initial estimates of those parameters. The 
program then calculates the least-squares parameter estimates. 
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the standard error of these estimates, the MA constant term, 
the residual variance, the autocorrelations and partial 
autocorrelations of the residuals, the final model sum 
of squared residuals and a Chi-square goodness-of-f it sta- 
tistic (the portmanteau test) . Additionally, it plots 
offline the autocorrelations and partial autocorrelations of 
the residuals. WMARQRDT uses the IMSL subroutines FTRDIF, 
FTAUTO, LINVIF and MDCDFI, the SSPLIB routine PLOTS, and 
the Time Series Editor resident subroutines PARSH, MARQRT, 
SUMSQ, SVJAPB and FORMB. The calculations are performed in 
double-precision arithmetic. 

J. XSUMSQ PROGRAM 

The XSUMSQ program accpets an already transformed and 
differenced series in FILE FT02F001 and general seasonal 
and non-seaonal Box-Jenkins parameter values for a general 
seasonal model from the terminal. The program then calcu- 
lates the residual sum of squares for these parameters, 
giving the user a feel for the "goodness-of-f it" the 
parameter estimates provide in modeling the given time 
series. XSUMSQ uses the Time Series Editor resident sub- 
program XSUMSQ. All output is printed at the terminal. 

The calculations are performed in double-precision arithmetic. 

K. FORECAST PROGRAM 

The FORECAST program uses the estimated Box-Jenkins 
(seasonal or non-seasonal) model to compute forecasts of 
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the transformed and differenced time series. FORECAST re- 
quires that the time series to be forecasted reside in FILE 
FT02F001 in transformed form. During execution the input 
required at the terminal includes the forecast origin, the 
numbers and estimated values of the AR and MA parameters 
of the model, the overall MA constant, the plot origin 
index, the maximum forecast lead time, the order of differ- 
encing in the model, and a level of significance for the 
forecast confidence limits. The program needs the transfor- 
mation parameters from FILE FT07F001, created when using 
program TRANS, and performs the inverse transform to return 
the data to its original form. FORECAST also undiffer- 
ences differenced data. The program output includes the 
forecasts up to the given maximum lead time, the deviations 
from each forecast for the confidence limits, and plots of 
forecasts and confidence limits plotted offline. The 
FORECAST program uses IMSL subroutine FTCAST and the special IBM 
SSPLIB subroutine UTPLT8. The plot itself and output data 
are written onto FILE FT08F001 on the user's disc; since 
this plot file takes considerable space, it is usually 
best to erase it after execution of the program is completed 
and the plot itself has been printed offline. 

L. ROOTS PROGRAM 

The ROOTS program is used to determine the roots of 
the characteristic equation for a general ARIMA model. The 
roots are useful for testing for stationarity and for 
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determining the form of the forecast function. ROOTS uses 
the IMSL subroutine ZPOLR to calculate the roots; it requires 
as input the number of AR parameters in undifferenced form,^ 
and the actual values of these parameters, as determined 
from the ESTIMATE or WMARQRDT programs. The roots are 
printed out at the terminal. If desired, roots can also 
be obtained for the MA polynomial, simply by substituting 
the appropriate MA values for the AR values as input. 

M. GENERATE AND SIMULATE PROGRAMS 

Because of their similarities, the GENERATE and SIMULATE 
programs are described together. The GENERATE program per- 
mits the user to generate a time series from any non- seasonal 
ARIMA model he specifies. The user must identify the model 
and give values for the model parameters and starting con- 
ditions; a random number seed must also be input. The 
program takes the specified model, generates random noise 
terms, and calculates as many values of the time series as 
desired. The GENERATE program can be useful for purposes 
of classroom instruction in the generation of a wide variety 
of time series examples for model identification exercises. 



^Suppose the model was ARIMA (1,1,0). The differenced form 
would be (1-(J), B) V^z . =9 + a.. Then, the undifferenced form 

would be written (1-(1+4),)B + <t),B^)z^ = 0 + a, , where the 

J. X O u 

left-hand side was found by multiplying (l-<t>^B) by (1-B) = V . 
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It could also be useful in the diagnostic phase of model 
checkout; a time series could be generated from the esti- 
mated model, and its properties compared with those of 
the original series. If large discrepancies occur in this 
comparison,' this may be evidence that the model is inade- 
quate. The model output, the generated time series, is 
written in FILE FT02F001 as well as being printed offline. 

The GENERATE program uses IMSL subroutine FTGENl. 

The SIMULATE program provides the capability of gen- 
erating any number of simulated time series. This program 
is useful for prediction of what might happen in the future; 
it also demonstrates that within a given model, the time 
series actually observed can vary considerably. This program 
uses GENERATE, but also requires as input the actual series 
in FILE FT02F001^ and the number of simulated series the 
user wishes to generate. The series should be already 
transformed and differenced. Additionally, SIMULATE allows 
the user to select values of the original time series as 
starting values for the simulated series. The output con- 
sists of the simulated series printed at the terminal. The 
SIMULATE program uses IMSL subroutine FTGENl. 

N. HELP OPTION 

The HELP option is a modification to the Time Series 
Editor CP/CMS program itself that allows the user to obtain 



^If desired, the user may choose to input his own starting 
values, and not use the actual time series values. 
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information about Editor programs after the initial 
introductory phase has been completed. 

0. SUMMARY 

This chapter has provided brief descriptions of the 
options available for analysis of time series data in the 
Time Series Editor. The following two chapters contain 
detailed examples of time series analysis for both a non- 
seasonal and a seasonal series; the accompanying terminal 
sessions and computer output for analysis of these series 
are contained in full in Appendices B and C. Appendix A 
to this report contains a User's Guide for the Time Series 
Editor, while Appendix D contains a complete listing of 
the Time Series Editor and resident programs. Appendix E 
contains listings of the data sets used in the analyses of 
the next two chapters. 
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V. EXAMPLE NON-SEASONAL TIME SERIES ANALYSIS 



This chapter provides a description of a complete analy- 
sis of a non-seasonal time series using the Time Series 
Editor. The time series analyzed is Series C (Chemical 
Process Temperature Readings), taken from Box and Jenkins 
[Ref. 4, p. 528]. This time series was selected since it 
is fully analyzed and discussed in Reference 4, so that the 
interested user can compare for himself the results pre- 
sented there with the results given by the Time Series 
Editor. After reading this overview of the analysis, the 
reader can follow the actual user session and associated 
user session Editor program output contained in Appendix B 
to this report. 

The first step in the analysis of series C is to plot 
the series using program PLOT. The plot reveals rather wide 
fluctuations in the series, but not the sort of explosive 
nonstationary behavior that would render an attempt to model 
fruitless. The plot also reveals that the time series has 
a large amount of momentum (movements of the series tend to 
resist changes of direction over time) . This is character- 
istic of ARIMA models with two or more differences. 

As a second step, the autocorrelations (autos) and 
partial autocorrelations (pautos) , mean, and variance of 
the series were estimated using AUTO. The plot shows that 
the autos dampen out slowly in a nearly linear fashion. 
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This is an indication that the process is nonstationary 

and that one or more differences are required to make it 

stationary. The pauto plot is not informative when the 

autos fail to dampen out rapidly. 

As suggested by the plots of the original series and 

its auto and pautos, the program DIFF was used to examine 

the series with both one and two non-seasonal differences; 

1 2 

i.e. , the series = V z. and z" = V z. were derived. 

t t t t 

The plots suggest that two differences might be required to 
achieve stationarity in the series. However, one might be 
able to get by with a single difference and an AR parameter 
near unity, since the autos of z^ dampen out slowly. Thus, 
two possible model candidates are suggested; 

1. ARIMA(1,1,0) : (1 - 4>^B)V^z^ = , 

and 

2. ARIMA(0,2,2) : V^z^ = 9^ + (1 - Q ' 

For purposes of estimation, the second model was extended 
to two moving average parameters. Such overfitting is often 
done to see if the estimated, moving average parameters turn 
out to be near zero, thus confirming the tentative identification. 

The next step in the analysis process is to calculate 
the maximum likelihood estimates of the model parameters. 
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The program ESTIMATE was used to do this. The estimated 
parameters for the ARIMA( 1,1,0) model were; 

= 0.8073 

0Q = -0.006789 

a.^ = 0.0177633 

a 

The autos and pautos of the residuals were calculated 

to test the model . The correlations should appear to be 

estimates of a pure white noise process if the model is 

adequate; the (1,1,0) model seemed to pass that test. The 

2 

chi-square goodness-of-f it test results in a x =28.88 with 
24 df and a signficance level of 0.2247. Thus, there is 
no strong evidence tc suggest that the (1,1,0) model is 
inadequate. The model can be written, in undifferenced 
form: 

(1 - 1.8073B + 0.8073B^)z^ = a^ . 

The 6^ is omitted here since it is nearly zero in value. 

The estimated parameters for the (0,2,2) model were 

= 0.1378 

02 = 0.1296 

0 = -0.0026787 

o 

0 ^ = 0.0189559 

As before, the correlation plots of the residuals fail to 

suggest any inadequacy of the model. However, the chi-square 

2 

goodness-of-f it test for this model yielded ax =36.75 
with 23 df and a significance level of 0.0346. Therefore, 
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if the ARIMA(0,2,2) model were a "correct" one, a chi-square 
value as large as 36.75 would occur by chance with a proba- 
bility of only 0.0346. Due to its simplicity (parsimony 
is a desirable feature in time series models) and better 
fit, the ARIMA(1,1,0) model was chosen as the "best" of 
the two alternatives originally selected. 

The ARIMA( 1,1,0) model calculated was then used to make 
forecasts and confidence limits for these forecasts, using 
program FORECAST. The forecasts were made for 25 periods 
into the future, with plot origin 200 and forecast origin 
220. The length of the series C itself is 226; with one 
difference, it becomes 225. The plot in Appendix B shows 
the forecasted values, along with the confidence limits at 
level of significance 0.10; the plot shows that the width 
of the confidence limits increases substantially as the lag 
gets large. 

Finally the program ROOTS was used to calculate the 
roots of the characteristic equation for the ARIMA( 1,1,0) 
model, which can be written: 

1 - 1.8073B - 0.8073B^ = 0 

The roots calculated by ROOTS were 1.239 and 1.0. Also, 
the program SIMULATE was used to generate a simulated series 
from the ARIMA( 1,1,0) model, using the last two values of 
the actual series as starting values. The entire session 
lasted about two hours, including checkout of plots and 
consumed less than a minute of CPU time. 
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VI. EXAMPLE SEASONAL TIME SERIES ANALYSIS 



This chapter provides a description of a complete analy- 
sis of a seasonal time series using the Box-Jenkins tech- 
nique through the Time Series Editor. The time series 
analyzed is Series G (International Airline Passengers: 
Monthly Totals (Thousands of Passengers) January 1949 - 
December 1960), taken from Box and Jenkins [Ref. 4, p. 531]. 
This time series was selected since it was fully analyzed 
and discussed in Chapter 9 of Reference 4; again, the 
interested reader can compare for himself the results pre- 
sented by Box and Jenkins with those given by the Time 
Series Editor programs. After reading this overview of the 
analysis, the reader can follow the actual user session and 
associated Editor program output contained in Appendix C 
to this report. 

The plot shows both annual seasonality and an increasing 
trend. This suggests that a seasonal model with differencing 
may be needed. The autocorrelations (autos) and partial 
autocorrelations (pautos) were then calculated and plotted 
using AUTO. The plot of the autos shows very slow damping 
out with peaks at intervals of 12, indicating strong non- 
stationarity and a seasonal period of 12 months. Again, 
here the pauto plot is not informative, due to the slow 
damping out of the autos. As it is sometimes useful in 
dealing with seasonal models, a natural log transform was 
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made on the series G data using program TRANS. The plot 
of the logged data shows little change in non-stationarity , 
and the autos and pautos of the logged data still indicate 
seasonality and non-stationarity. 

Next, differencing to achieve stationer ity was tried, 
using program DIFF; both one seasonal and one non-seasonal 
difference of the logged data was taken. The autos and 
pautos of the data '^^^£2 2 ^) appear much improved, 

rapidly damping out. They also exhibit sharp peaks at the 
periods of seasonality, indicating the need for seasonal 
parameterization in the model. Since positive identification 
of the model type from the auto and pauto plots was not 

possible, two were postulated as close candidates: 

« 

1. ARIMA (0,1,1) X (0,1,1)^^ , written: 

V^V^'-^dn z^) = + (1-0^B) ' 



and 



2. ARIMA (1,1,1) X (0,1,1)^^ , written: 

(1-9^B) V^V£'-2(ln z^) = 9^ + (1-9 j_B) ( 1-0^ ^ ^t * 

The next step in the analysis was to calculate the 
parameter estimates for both models. YESTSEAS was used to 
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obtain initial parameter estimates. For the model 

12 

ARIMA (0,1,1) X (0,1/1) , YESTSEAS yielded starting values 



6, , , = 0.390425 

1 (s) 



®l,12(s) 



0.534397 . 



For this model, calculation of final parameter estimates 
in WMARQRDT using the starting values above yielded the 
following in seven iterations: 



= 0.377152 



®1 12 ^ 0.572387 



6 = 0.000291 ; 

o 



of fit of the 
a probability 

YESTSEAS calcu- 

‘^l(s) ^ 0.112699 

9^(s) = 0.490883 



the chi-square statistic for residual lack 
2 

model was x ~ 29.713571, with df = 38 and 

2 

of exceeding the x value of .829. 

For the ARIMA (1,1,1) x (0,1,1)^^ model 
lated starting values: 
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0.533679 ; 



®l,12(s) 

for this model, WMARQRDT yielded the following parameter 
values in ten iterations: 

= 0.146667 

0^ = 0.509816 

©1 ^2 ^ 0.573388 . 

The chi-square statistic for residual lack of fit of the 

model was ~ 28.302543, with df = 37 and a probability 

2 

of exceeding the x value of .847. 

Next, the plots of autos and pautos of the residual for 
both models were examined; both models have autos and pautos 
that dampen out rapidly, and essentially appear to be 
"white noise", with the exception of a small peak at auto 
and pauto of lag 12, corresponding to the seasonal period. 
The chi-square statistics indicate that both models have 
"about" the same goodness of fit to the data. Consequently, 
both models were used to make forecasts of the series. 

Program FORECAST was used to forecast series G into the 
future for 30 time periods, with a forecast origin of 131; 
the plot origin selected was index 100. As expected, both 
models do about the same job of forecasting the series. 
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The significance level for the confidence limits for both 
plots was 90%. 

As a final check of both models, the program XSUMSQ 
was used to calculate the sum of squared residuals for 
models. The results were as follows: 

1. For ARIMA (0,1,1) x (0,1,1)^^ , 

S^(0„, 0. A = 0.181929 

O J. X / XZ a. 

2. For ARIMA (1,1,1) X (0,1,1) , 

0o'®l'®l,12'^a^ "" 0.181643 . 

It can be seen here that both models are nearly equally 
"good", and either could be used with a reasonable amount 
of confidence for forecasting series G. 
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VII. SUMMARY 



There has been an increasing need in military applica- . 
tionS/ as well as in the general business area, for a 
computer-aided time serie? analysis and forecasting capability 
that is both powerful and easy to use. Current military- 
associated functions in which time series techniques are 
valuable include material inventory management, recruiting, 
personnel management, budget analysis, force level projections, 
and perhaps even short-term weather prediction. 

The Time Series Editor that has been described in this 
report provides a unified collection of programs that greatly 
facilitate the conduct of a complete time series analysis 
and forecasting evolution using Box-Jenkins methodology. 

The routines guide the user from the data input stage, 
through the data analysis and parameter estimation stages, 
and finally through the model building, diagnostic checking 
and forecasting stages. With its simple input requirements, 
the Editor can be quickly mastered by even a beginning 
computer user, who has a basic understanding of Box-Jenkins 
methodology. 

The Box-Jenkins methodology has been described broadly 
in Chapter II and descriptions of the programs resident 
in the Time Series Editor have been given in Chapter IV. 

Example analyses of both non-seasonal and seasonal time 
series are provided in Chapters V and VI, with complete 
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user sessions and output packages for these analyses 
included as Appendices B and C. The Appendices also con- 
tain a Time Series Editor User's Guide, as well as complete 
program listings. 

Although the Time Series Editor in its present form 
provides an excellent capability for interactive time 
series analysis using the Box- Jenkins methodology, there 
are additions to the program set that might be made in the 
future which would improve the Editor's capability and 
utility. It is first recommended that the model diagnostic 
testing capability be extended to include a periodogram 
analysis, and/or other tests related to the spectral analy- 
sis of time series [see Ref. 4, p. 294]. Another useful 
addition might be to include in the Time Series Editor a 
package capable of improving the ease of entering data 
and simplifying the input error correction problem. 
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APPENDIX A 



USER'S GUIDE TO THE TIME SERIES EDITOR 

The Time Series Editor is a collection of FORTRAN 
programs driven by a control program called TIMESER EXEC 
in the CP/CMS executive language that has been designed 
specifically for the analysis and forecasting of time 
series data using the general Box- Jenkins methodology. 

This Guide contains information essential for the user to 
access the Time Series Editor, enter time series data, 
build a model and obtain the desired output. 

I. DATA INPUT 

The Time Series Editor requires that time series data 
be entered into the sequential FORTRAN input/output file 
named FILE FT02F001. This can be done either online or 
via cards read offline. The following diagram illustrates 
the proper card deck arrangement. When your deck is ready 
for input, give it to a computer center system operator 
to direct it to your disc space. 

If the user has a data deck already punched up in a 
format other than 5F15.6, he may enter the series as above 
into FILE FT03F001 ( without the length of series card) and 
use the Editor program ZFORMAT to transform it into the 
FILE FT02F001 in the proper format, without destroying his 
original file. This program is described in the next section. 
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DATA ENTRY DECK 




ro 
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II. TABLE OF OPTIONS 

This Table provides the user with the basic information 
necessary to understand the data requirements, functions 
and output options for each program in the Editor. 
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III. THE BASIC USER SESSION 



To use the Time Series Editor, the user must log into 
CMS, get into CP, link to the disc storage area where the 
Time Series Editor resides, reimplement CMS, log into the 
general user and Time Series Editor disc areas, and enter 
the TIMESER routine. This section will provide explicit 
guidelines to enable the user to perform the above steps 
on the NPS CP/CMS system. Commands marked with an asterisk 
(*) are those actually entered on the terminal by the user 
(the asterisk itself is omitted) . Those without an aster- 
isk and those written in all capital letters are system 
responses at the terminal. Numbered sentences are comments, 
which will not appear during an actual user session. The 
instructions and system responses assure the user is on an 
IBM 2741 Input/Output Terminal. Some minor modifications 
may be necessary if other terminals are used. 



1. Turn the terminal on, depress the BREAK key, and wait 
for the system to respond: 

CP-67 online xd.65 gsyosu 

2. Depress the ATTN key. The roll bar will advance and 
the keyboard will unlock. Then enter: 

*login aaaapbb 450k 

1 

3. aaaa is the user's identification number, and nn is the 
terminal niomber (usually written on the terminal) . 

For example, if the user's ID number is 1621 and the 
terminal number is 44, the input would be: 
login 1621p44 450k. The addition of 450 k to the normal 
login command is necessary to execute the program 
WMARQRDT in the Editor; for users not planning to 
execute this program during a session, this addition 
is not necessary. 
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4. The system will respond with the statement: 

ENTER PASSWORD : 

5. The user then enters his password, or the general users 
password npg; 

*pas sword 

6. The system will then respond; 

ENTER 4-DIGIT PROJECT NUMBER FOLLOWED BY 4-CHARACTER 
COST CENTER CODE: 

7. The user then enters: 

*gggghhhh 

8. gggg is the assigned project number, and hhhh is the 
user's section designator or the faculty code. 

9. The system will respond with the message of the day, 
such as: 

CP/CMS HOURS ... 0930=2200 (MON-THURS) ... 0930-1800 (FRI) 
•OUTPUT RETAINED 5 DAYS 
Cms Version 3.25 

10. At this point, the user is in CMS. He must then get into 
CP; this can be done by hitting the ATTN key. The system 
will then respond; 

CP 

11. The user must then link to the TIME SERIES EDITOR; 
this is accomplished by entering: 

*link 2069p 191 192 

12. The system will respond with: 

ENTER PASSWORD: 

13. The password (read only) to enter the Editor is: 

*timser 

14. The system then responds: 

SET TO READ ONLY 
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15. The user now implements CMS by: 

*ipl cms 

16. The system will respond: 

CMS Version 3.25 

17. Now the user must log into both the general user and 
the Time Series Editor area by entering: 

*login 191 

18. The system will respond with a message such as: 

R; 

19. The user then enters the command: 

♦login 192 t,p 

20. The system will respond: 

. T (192) R/0 

R; 

21. The user can then enter the Time Series Editor 
(guided version) by entering the command: 

♦timeser 

22. The system will respond: 

EACH 2 SECONDS EXECUTION TIME IS INDICATED BY * 

YOU HAVE ENTERED THE TIME SERIES EDITOR 

PLEASE RESPOND TO EACH QUERY WITH AN INPUT AT THE TERMINAL. 
ENTER ONLY THE FIRST LETTER FOR A WORD RESPONSE. 

ENTER NUMERICAL VALUES VIA FORTRAN FORMAT.. 

TYPE INTEGER VALUES (RIGHT JUSTIFIED) FOR NAMES STARTING 
WITH I THROUGH N. TYPE FLOATING VALUES WITH DECIMAL FOR 
ALL OTHERS. 

DO YOU WANT A LIST OF THE OPTIONS? 

23. The user is then on his own, guided by the Exec routine. 

See the notes that appear at the end of this guide for 
additional information. Eventually the user will be 
asked: 

DO YOU WANT TO TRY AGAIN? 
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If a yes response is given, another sequence will begin; 
if the response is no, the user will be taken out of 
the Time Series Editor environment and returned to CMS. 
The system response will be: 

CONTROL RETURNED TO CMS 
R; 



25. The user can then log out of CMS by typing; 

*cp logout 

26. The system will respond with; 

CONNECT= 00:08:02 VIRTCPU= 000:07.98 TOTCPU= 000.10.94 
LOGOUT AT 14.22.04 on 10/16/78 

27. The user should then turn off his terminal and tear off 
the output from his session. 

The more experienced user can dispense with the 
"welcome aboard" section of the Time Series Editor and get 
right down to business by using the shortened version of the 
Editor. This shortened version may be entered by linking in 
the normal way, and then entering the Editor by typing the 
COMMAND 



*timeser s (asterisk omitted) . 



The system will immediately respond: 



ENTER LETTER FOR OPTION YOU WANT. 



The session inside the Editor then begins. 
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IV. BRIEF SAMPLE USER SESSION 



A brief sample user session is given below; it includes 
copies of the offline output generated during the session. 
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repeat login nur@pn§ 



login 1621p44 450k 
ENTER PASSWORD; 

ENTER 4-DIGIT, PROJECT NUMBER FOLLOWED BY 4-CHARACTER COST CENTER CODE 
0444tl72 * 

READY AT 17.24.38 ON 09/16/78 
CMS Version 3.25 

stat 

P (191): 29 FILES; 241 REC IN USE, 55 LEFT (of 296), 81% FULL (2 CYL) 
R; 

cp q f 

FILES:- NO RDR, NO PRT, NO PUN 
R; 

CP 

link 2069p 191 192 
ENTER PASSWORD: 

SET TO READ ONLY 

ipl cms 

CMS Version 3.25 

login 191 

R; 

login 192 t,p 
T (192) R/0 
R; 

timeser 

EACH 2 SECONDS EXECUTION TIME IS INDICATED BY *. 

YOU HAVE ENTERED THE TIME SERIES EDITOR. 

PLEASE RESPOND TO EACH QUERY WITH AN INPUT AT THE TERMINAL. 

ENTER ONLY THE FIRST LETTER FOR A WORD RESPONSE. 

ENTER NUMERICAL VALUES VIA FORTRAN FORMAT. 

TYPE INTEGER VALUES (RIGHT JUSTIFIED) FOR NAMES STARTING 

WITH I THRU N. TYPE FLOATING VALUES WITH DECIMAL FOR ALL OTHERS. 

DO YOU WANT A LIST OF THE OPTIONS? 

y 



OPTION DESCRIPTION 

GENERATE GENERATE ANY ARIMA TIME SERIES 

AUTO CALCULATE AUTOCORRELATIONS, PAUTOS, MEAN AND VARIANCE 

PLOT PLOT A TIME SERIES 

ESTIMATE CALCULATE MAX LIKELIHOOD ESTIMATES OF ARMA PARAMETERS 

DIFF DIFFERENCE A TIME SERIES 

FORECAST FORECAST FUTURE VALUES, CONSTRUCT CONFIDENCE INTERVALS 

TRANS -“TRANSFORMS VALUES OF A TIME SERIES 

ROOTS DETERMINES ROOTS OF ARIMA CHARACTERISTIC EQUATION 

ZFORMAT ALTER DATA FILE TO FORMAT 5F15.6 

CMSWORK PERFORM CP/CMS COMMANDS IN TIMESER EXEC 

SIMULATE SIMULATE NONSEASONAL TIME SERIES 

YESTSEAS CALCULATE INITIAL SEASONAL PARAMETERS 

WMARQRDT HARQUARDT SOLUTION FOR PARAMETER ESTIMATES 

XSUMSQ CALCULATE SUM OF SQUARES FOR ARBITRARY PARAMETERS 



WOULD YOU LIKE MORE INFO? 

y 
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ENTER OPTION YOU WANT INFO ABOUT, 
a 

AUTO THIS PROGRAM CALCULATES AUTOCORRELATIONS, PARTIAL 

AUTOCORRELATIONS, THE MEAN AND THE VARIANCE FOR A GIVEN TIME SERIES 
WHICH MUST RESIDE IN FILE FT02F001. THE PROGRAM USES 
FTAUTO IN THE IMSL LIBRARY. THE AUTOCORRELATIONS AND PAUTOS CAN BE 
PLOTTED OFFLINE. 



DO YOU WANT INFO ABOUT ANOTHER OPTION? 
n 

DO YOU WANT TO TRY A SESSION? 

y 

ENTER LETTER FOR OPTION YOU WANT. 

C 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 



listf * 


ft02f001 








FILENAME 


FILETYPB 


MODE 


NO. REC. 


DATE 


SERG 


FT02F001 


PI 


3 


9/16 


SERC 


FT02F001 


PI 


5 


9/16 


LNSERG 


FT02F001 


PI 


3 


9/16 


FILE 


FT02F001 


PI 


1 


9/16 



erase file ft02f001 

alter sere ft02f001 ?1 file ft02f001 pi 
stat 

P (191); 28 FILES; 240 REC IN USE, 56 LEFT (of 296), 81% FULL (2 CYL) 
&goto -ques 

DO YOU WANT TO GO AGAIN? 

y 



ENTER LETTER FOR OPTION YOU WANT, 
a 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

AUTOCORRELATIONS 

0.978 0.944 0.902 0.854 0.802 0.748 0.692 0.635 0.579 0.323 

0.468 0.413 0.359 0.305 0.253 0.201 0.150 0.098 0.047 -0.003 

-0.052 -0.101 -0.151 -0.200 -0.248 



PARTIAL AUTOCORRELATIONS 

0.978 -0.260 -0.157 -0.093 -0.058 -0.045 -0.012 -0.038 -0.022 -0.010 
-0.036 -0.041 -0.038 -0.024 -0.037 -0.027 -0.032 -0.070 -0.048 -0.024 
-0.034 -0.061 -0.079 -0.048 -0.037 



MEAN= 22.9739 VARIANCE * 4.22273 

ENTER TITLE FOR PLOTS, 
autos and pautos for series c data 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT. 

C 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 

WHEN FINISHED TYPE: &GOTO -QUES 
offline print file ft02f001 
alter file ft02f001 pi sere ft02f001 pi 
&goto -ques 

DO YOU WANT TO GO AGAIN? 
n 

CONTROL RETURNED TO CMS 
R? 
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AUTCCCRBFLATICNS with 2 SIGHA BANDS 
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V. PROBLEM CONTROL NOTES 



This section will cover corrective measures that can be 
taken when things fail to go as expected while in the 
TIMESER environment. 

a. A typing error in the CMS environment can be 

corrected by typing the § character as many times as is 
required to back up and then type the correct values. For 
example, if the user typed timesre , the user could correct 
the mistake by typing two @ signs, followed by the correct 
spelling er, as follows: timesre@@re . An entire line can 

be deleted by typing the character (or [ on some terminals) 

b. When working with TIMESER executive programs, the 
user should exercise care before hitting the return key (or 
control s on some terminals). If an input value is required 
and the return key is hit before the proper response is 
entered, the user will be likely to get thrown out of the 
editor and have to begin again. In most cases, errors can 
be corrected only before the return key is struck (in some 
programs you get a second chance for input) . 

c. Particular care should be taken for integer value 
input, which must be right justified in the format field. 

The editor will advise the user in all cases where the 
integer format is other than II. 

d. If for any reason the user finds himself in a debug 
or error condition (caused by erroneous data, or a "blowup" 
in one of the non-linear optimization routines usually 
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caused by very poor initial input values) , the following 
procedure will get the user back into the normal TIMESER 
environment; 

(1) depress the ATTN key twice; this gets 
the user into CP; hit the ATTN key again, 
and then type kx to kill the execution. 

(2) re-ipl CMS, and login 191 and then login 
192 t,p. 

(3) then the user can type TIMESER or TIMESER S, 
and return to the TIMESER environment 

On the next page is a sample user session where an error 

causing a debug condition has occurred. 

Here the user executed a program that required data 
in FILE FT02F001, and the FILE did not exist. As the 
example shows, recovery is quick. Simply hit the break 
button to get into CP, login 191, login 192 t,p, and then 
type timeser s to return immediately to the Editor environment. 
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ENTER LETTER FOR OPTION YOU WANT 
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APPENDIX B 



KEYBOARD SESSION AND OUTPUT FROM ANALYSIS OF 
TIME SERIES C (NON-SEASONAL) 

This appendix contains the complete terminal keyboard 
session and all offline output from the analysis of Box 
and Jenkins' [Ref. 4] time series C (Chemical process 
temperature readings) using the Time Series Editor. This 
appendix is intended to supplement the discussion of the 
analysis of series C discussed in Chapter V of this report. 

The terminal listings are presented first, followed by 
the computer offline output. 
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repeat login nur0pn@ 



login 1621p44 400k 
ENTER PASSWORD; 

ENTER 4-DIGIT PROJECT NUMBER FOLLOWED BY 4-CHARACTER COST CENTER CODE 
0444tl72 

READY AT 14.23.02 ON 09/16/78 
CMS Vetsiorf 3.25 



listf * 
FILENAME 


ft02f001 

FILETYPE 


MODE 


NO.REC. 


DATE 


FILE 


FT02F001 


PI 


3 


9/15 


SERG 


FT02F001 


PI 


3 


9/16 


SERC 


PT02F001 


Pi 


5 


9/16 



R; 

erase file £t02£001 

R* 

alter setc £t02£001 pi file £t02f001 pi 
R; 

CP 

link 2069p 191 192 
ENTER PASSWORD: 

SET TO READ ONLY’ 

ipl cms 

CMS Version 3.25 

login 191 
R; 

login 192 t,p 
T (192) R/0 
R; 



..er 

-aCH 2 SECONDS EXECUTION TIME IS INDICATED BY *. 

YOU HAVE ENTERED THE TIME SERIES EDITOR. 

PLEASE RESPOND TO EACH QUERY WITH AN INPUT AT THE TERMINAL. 

ENTER ONLY THE FIRST LETTER FOR A WORD RESPONSE. 

ENTER NUMERICAL VALUES VIA FORTRAN FORMAT. 

TYPE INTEGER VALUES (RIGHT JUSTIFIED) FOR NAMES STARTING 

WITH I THRU N. TYPE FLOATING VALUES WITH DECIMAL FOR ALL OTHERS. 

DO YOU WANT A LIST OF THE OPTIONS? 

y 



OPTION DESCRIPTION 

generate -GENERATE ANY ARIMA TIME SERIES 

AUTO CALCULATE AUTOCORRELATIONS, PAUTOS, MEAN AND VARIANCE 

PLOT PLOT A TIME SERIES 

ESTIMATE CALCULATE MAX LIKELIHOOD ESTIMATES OF ARMA PARAMETERS 

DIFF DIFFERENCE A TIME SERIES 

FORECAST FORECAST FUTURE VALUES, CONSTRUCT CONFIDENCE INTERVALS 

TRANS TRANSFORMS VALUES OF A TIME SERIES 

ROOTS DETERMINES ROOTS OF ARIMA CHARACTERISTIC EQUATION 

ZFORMAT ALTER DATA FILE TO FORMAT 5F15.6 

CMSWORK PERFORM CP/CMS COMMANDS IN TIMESER EXEC 

SIMULATE SIMULATE NONSEASONAL TIME SERIES 

YESTSEAS CALCULATE INITIAL SEASONAL PARAMETERS 

WMARQRDT MARQUARDT SOLUTION FOR PARAMETER ESTIMATES 

XSUMSQ CALCULATE SUM OF SQUARES FOR ARBITRARY PARAMETERS 



WOULD YOU LIKE MORE INFO? 

y 
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ENTER OPTION YOU WANT INFO ABOUT, 
a 

AUTO THIS PROGRAM CALCULATES AUTOCORRELATIONS, PARTIAL 

AUTOCORRELATIONS, THE MEAN AND THE VARIANCE FOR A GIVEN TIME SERIES 
WHICH MUST RESIDE IN FILE FT02F001. THE PROGRAM USES 
FTAUTO IN THE IMSL LIBRARY. THE AUTOCORRELATIONS AND PAUTOS CAN BE 
PLOTTED OFFLINE. 



DO YOU WANT INFO ABOUT ANOTHER OPTION? 

y 

ENTER OPTION YOU WANT INFO ABOUT, 
d 

DIFF THIS PROGRAM TAKES A SERIES IN FILE FT02F001, ASKS IF 

IT IS SEASONAL, AND THEN PERFORMS NONSEASONAL AND/OR SEASONAL 
DIFFERENCING, AS REQUESTED. THE ORIGINAL DATA REMAINS IN FILE 
FT02F001, WHILE THE DIFFERENCED SERIES IS PUT IN PILE FT03F001. 



DO YOU WANT INFO ABOUT ANOTHER OPTION? 
n 

DO YOU WANT TO TRY A SESSION? 

y 

ENTER LETTER FOR OPTION YOU WANT, 
a 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

AUTOCORRELATIONS 



0.978 0.944 

0.468 0.413 

-0.052 -0.101 



0.902 0.854 

0.359 0.305 

-0.151 -0.200 



0.802 0.748 

0.253 0.201 

-0.248 



0.692 0.635 

0.150 0.098 



0.579 0.523 

0.047 -0.003 



PARTIAL AUTOCORRELATIONS 

0.978 -0.260 -0.157 -0.093 -0.058 
0.036 -0.041 -0.038 -0.024 -0.037 
0.034 -0.061 -0.079 -0.048 -0.037 



-0.045 -0.012 -0.038 -0.022 -0.030 
-0.027 -0.032 -0.070 -0.048 -0.024 



MEAN» 22.9739 VARIANCE = 4.22273 

ENTER TITLE FOR PLOTS. 

autos and pautos of series c data / undifferenced 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

♦DO YOU WANT TO GO AGAIN? 

y 
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ENTER LETTER FOR OPTION YOU WANT, 
d • ‘ 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS.,. 

IS YOUR TIME SERIES SEASONAL? 
n 

ENTER NUMBER OF NONSEASONAL DIFFERENCES. 

1 

♦DO YOU WANT TO PLOT AUTO AND PAUTO OP TRANSFORMED DATA? 

y 

EXECUTION BEGINS... 

AUTOCORRELATIONS 



0.805 

0.097 

0.044 


0.653 

0.094 

0.048 


0,526 

0,074 

0.000 


0.442 0.380 

0,073 0,070 

-0.052 -0.097 


0,318 

0.072 


0.262 

0.089 


0.186 

0.048 


0.139 

0.041 


0.144 

0.040 




PARTIAL AUTOCORRELATIONS 












0.805 

•0.137 

0.038 


0.010 

0.094 

-0.022 


-0.007 

-0.027 

-0.116 


0.051 0,027 

0.034 0.006 

-0,076 -0.008 


-0.019 

0.012 


-0.013 

0,043 


-0,080 

-0.121 


0,020 

0.060 


0,117 

-0.005 



MEAN=-,346667E-01 VARIANCE » 0.531982E-01 
ENTER TITLE FOR PLOTS. 

autos and pautos of series c data / one nonseasonal difference 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 



DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT, 
d 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

IS YOUR TIME SERIES SEASONAL? 
n 

ENTER NUMBER OF NONSEASONAL DIFFERENCES. 

2 

DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 

y 

EXECUTION BEGINS,, . 

AUTOCORRELATIONS 

-0.079 -0.065 -0.122 -0.063 0.013 -0.018 0.049 -0.052 -0.124 0.122 
-0.122 0,072 -0,077 0,029 -0.011 -0.058 0.171 -0.101 -0,013 -0,020 
-0.007 0,135 0,014 -0,013 -0,135 



PARTIAL AUTOCORRELATIONS 



-0.079 

-0.143 

0.004 



-0.072 

0,023 

0.081 



-0.135 

-0.089 

0.056 



-0.094 

-0.005 

-0.041 



- 0.022 

-0.039 

-0,099 



-0.051 

-0.074 



0.022 -0.060 -0.145 0.094 

0.130 -0.105 -0.014 -0.053 



MEAN*-,267866E-02 VARIANCE = 0.198143E-01 



ENTER TITLE FOR PLOTS. 

autos and pautos of series c data / two nonseasonal differences 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 
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DO YOU WANT TO GO AGAIN? 

y 

* 

ENTER LETTER FOR OPTION YOU WANT* 

C 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE* 

WHEN FINISHED TYPE:, &GOTO -QUES 
alter file ft02f001 pi sere ft02f001 pi 
alter file ftOBfOOl pi file ft02f001 pi 
&gpto -*ques 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT, 
e 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS..* 

ENTER NUMBER OF AR PARAMETERS (EXCLUDE DIFFERENCES). 

0 

ENTER NUMBER OF MA PARAMETERS. 

2 

LENGTH OF TIME SERIES - 224 
0 AR PARAMETERS 2 MA PARAMETERS 



MA PARAMETERS ARE: 

THETA(1)» 0.1378 

THETA(2)== 0.1296 

MA CONSTANT* -0.26787E-02 WHITE NOISE VARIANCE* 0.189559E-01 



CHI SQUARE RESIDUAL LACK OF FIT VALUE * 36.75 DF » 23 

SIGNIFICANCE * 0.0346 

DO YOU WANT TO PLOT AUTO AND PAUTO OF RESIDUALS? 

y 

♦EXECUTION BEGINS... 

AUTOCORRELATIONS 



0.019 0.034 

-0.142 0.069 

0.009 0.133 



-0.134 -0.093 
-0.098 0.025 

0.005 -0.015 



-0.012 -0.043 
0.004 -0.055 
-0.133 



0.005 -0.063 
0.157 -0.090 



-0.125 0.088 

0.005 -0.006 



PARTIAL AUTOCORRELATIONS 



0.019 0.033 

-0.166 0.017 

0.024 0.069 



-0.135 -0.090 
-0.104 -0.013 
0.024 -0.069 



0.001 -0.055 
-0.025 -0.089 
- 0.102 



-0.019 -0.071 
0.123 -0.123 



-0.143 0.087 

-0.021 -0.025 



MEAN»0.236684E-02 VARIANCE * 0.189503E-01 
ENTER TITLE FOR PLOTS. 

autos and pautos of residuals / arima(0,2,2) model of series c data 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

DO YOU WANT TO GO AGAIN? 

y 
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ENTER LETTER FOR OPTION YOU WANT, 
c 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 

WHEN FINISHED TYPE: SGOTO -QUES 

erase file £t02f001 

scat 

P (191): 27 PILES; 265 REC IN USE, 31 LEFT (of 296), 90% PULL (2 CYL) 
erase file* ftOSfOOl 

alter sere ft02f001 pi file ft02f001 pi 
sgoto -ques 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT, 
d 

*IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

IS YOUR TIME SERIES SEASONAL? 
n 

ENTER NUMBER OF NONSEASONAL DIFFERENCES. 

1 

DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 
n 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT. 

C 

‘ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 

WHEN FINISHED TYPE; &GOTO -QUES 
alter file ft02f001 pi sere ft02f001 pi 
alter file fc03f001 pi file ft02f001 pi 
&goto -ques 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT, 
e 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

ENTER NUMBER OF AR PARAMETERS (EXCLUDE DIFFERENCES). 

1 

ENTER NUMBER OF MA PARAMETERS. 

0 

* LENGTH OF TIME SERIES » 225 
1 AR PARAMETERS 0 MA PARAMETERS 

AR PARAMETERS ARE: 

PHKD- 0.8073 

MA CONSTANT- -0.66789E-02 WHITE NOISE VARIANCE- 0.177633E-01 

CHI SQUARE RESIDUAL LACK OF FIT VALUE - 28.88 DF - 24 

SIGNIFICANCE » 0.2247 ' 

DO YOU WANT TO PLOT AUTO AND PAUTO OF RESIDUALS? 

y 

EXECUTION BEGINS... 

AUTOCORRELATIONS 

0.016 0.014 -0.049 -0.005 0.061 0.024 0.078 -0.018 -0.087 0.130 
-0.093 0.084 -0.055 0.043 0.007 -0.040 0.169 -0.079 -0.001 -0.010 
0.000 0.129 0.014 -0.013 -0.127 
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PARTIAL AUTOCORRELATIONS 

0.016 0.014 -0.050 -0.004 0.063 0.019 0.075 -0.015 -0.087 0.14C 
-0.104 0.070 -0.048 0.042 0.003 -0.030 0.158 -0.088 0.013 -0.025 
0.030 0.086 0.039 -0.062 -0.101 



MEAN— .252517E-02 VARIANCE =• 0.177569E-01 
ENTER TITLE FOR PLOTS. 

autos and pautos o£ residuals / arima(l,l,0) model of series c data 

YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

DO YOU WANT TO GO AGAIN? 

y 

* 

ENTER LETTER FOR OPTION YOU WANT, 

c 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 

WHEN FINISHED TYPE: SGOTO -QUES 

erase file ft02f001 

alter sere ft02f001 pi file ft02f001 pi 
sgoto -ques 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT, 

f 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

IS YOUR SERIES SEASONAL? 

n 



ENTER NONSEASONAL DATA AS REQUESTED; 

ENTER NUMBER OF NON-SEASONAL DIFFERENCES. 

1 

ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 

1 

ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 

0 



NOW INPUT YOUR PARAMETER ESTIMATES: 



ENTER MA CONSTANT TERM, THETAO. 

-.00667 

ENTER NON-SEASONAL AR PARAMETER PHI(l). 

0.8073 

ENTER MAXIMUM FORECAST LEAD TIME VIA FORMAT 12. 

25 

ENTER INDEX FOR PLOT ORIGIN VIA FORMAT 13. 

200 

ENTER INDEX FOR FORECAST ORIGIN VIA FORMAT 13. 

220 

ENTER SIGNIFICANCE LEVEL FOR CONFIDENCE INTERVALS. 

0.10 

WAS YOUR DATA TRANSFORMED IN THE TRANS PROGRAM? 

n 



DO YOU WANT BASIC OUTPUT AT THE TERMINAL? 

y 



125 



I 



I 



I 



s 






THE LAST 10 WVEC VALUES 



0.222000E 


02 


0.218000E 


02 


0.213000E 


02 


0.208000E 


02 


0.202000E 


0.197000E 


02 


0.193000E 


02 


0.191000E 


02 


0.190000E 


02 


0.186000E 


THE 25 FORECAST VALUES; 














0.203896E 


02 


0.200517E 


02 


0.197722E 


02 


0.195399E 


02 


0.193456E 


0.191821E 


02 


0.190434E 


02 


0.189248E 


02 


0.188223E 


02 


0.187329E 


0.186541E 


02 


0.185837E 


02 


0.185203E 


02 


0.184623E 


02 


0.184089E 


0.183591E 


02 


0.183121E 


02 


0.182676E 


02 


0.182249E 


02 


0.181338E 


0.181439E 


02 


0.181050B 


02 


0.180669E 


02 


0.180295E 


02 


0.179926E 


THE 25 UPPER 1 


FORECAST CONFIDENCE LIMITS: 








0.206094E 


02 


0.205056E 


02 


0.204779E 


02 


0.205033E 


02 


0.205662E 


0.206551E 


02 


0.207621B 


02 


0.208812E 


02 


0.210079E 


02 


0.211390E 


0.212721E 


02 


0.214056E 


02 


0.215381E 


02 


0.216687E 


02 


0.217969E 


0.219221E 


02 


0.220442E 


02 


0.221630E 


02 


0.222783E 


02 


0.223903E 


0.224989E 


02 


0.226043B 


02 


0.227064E 


02 


0.228055E 


02 


0.229016E 


ALPHA FOR 


THE 


CONFIDENCE LIMITS IS: 0 


• 100 








DO YOU WANT TO 


GO AGAIN? 















y 

ENTER LETTER FOR OPTION YOU WANT, 
r 

EXECUTION BEGINS... 

SOLVES FOR ROOTS OF CHARACTERISTIC EQUATION. 

ENTER NO. OF AR PARAMETERS (UNDIFF. FORM). 

2 

ENTER AUTOREGRESSIVE PARAMETERS: 

ENTER PHI(l). 

1.8073 

ENTER PHI (2). 

-0.8073 

THE 2 ROOTS OF CHARACTERISTIC EQ. AND NORMS ARE: 
1.239 0.0 1.239 

1.000 0.0 1.000 

DO YOU WANT TO GO AGAIN? 

y 



ENTER LETTER FOR OPTION YOU WANT, 
d 

♦EXECUTION BEGINS... 

IS YOUR TIME SERIES SEASONAL? 
n 

ENTER NUMBER OF NONSEASONAL DIFFERENCES. 

0 

DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 
n 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT, 
s 

♦EXECUTION BEGINS... 

ENTER NUMBER OF AR PARAMETERS (UNDIFFERENCED FORM) 

2 

ENTER ESTIMATED AR PARAMETER PHI(l). 

1.8073 

ENTER ESTIMATED AR PARAMETER PHI (2). 

-0.8073 



02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 

02 
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AR PARAMETfiRS ARE 1.8073 -0.8073 

ARE THESE OK? 

y 

ENTER NUMBER OP MA PARAMETERS 

0 

ENTER OVERALL MA CONSTANT 
-0.00667 

ENTER ESTIMATED WHITE NOISE VAR 
0.0177 

DO you WANT STARTING VALUES TO BE THE LAST VALUES 
OF THE ACTUAL SERIES? 

y 

START(1)» 19.000000 

START'{2)» 18.799988 

STARTI NG' VALUE (S) OK? 

y 

ENTER RANDOM NUMBER SEED (BETWEEN 0 AND 1). . 

0.12345 

ENTER NUMBER OF VALUES YOU WANT TO SIMULATE VIA 13. 
50 

DO YOU WANT A PRINTOUT OF THE FIRST 20 VALUES? 

y 

. SIMULATED VALUES ARE: 



18.50224 


18.33957 


18.33694 


18.16435 


17.98804 


17.85420 


17.69398 


17.58267 


17.50722 


17.26422 


16.94655 


16.68689 


16.43227 


16;42499 


16.26595 


16.31497 


16.64629 


17.04015 


17.47137 


17.62828 


DO YOU WANT TO GO 


AGAIN? 








n 

CONTROL RETURNED 


TO CMS 









R; 
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PARTIAL ALTOCCBRELATIONS WITH 2 SIGMA BANDS. 
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SIGNIFICANCE LEVEL FOR CONFIOENCE IMERVALS: SO.OOO 

FCRLCAS7 CPIGINJ220 

MAXINLM FCFECAST LEAD TI^^E: 25 



APPENDIX C 



KEYBOARD SESSION AND OUTPUT FROM THE ANALYSIS OF 
TIME SERIES G (SEASONAL) 

This appendix contains the complete terminal keyboard 
session and all offline output from the analysis of Box 
and Jenkins' [Ref. 4] time series G (International airline 
passengers, monthly totals, Jan 1949 - Dec 1960) using the 
Time Series Editor. This appendix is intended to supplement 
the discussion of the analysis of series G given in Chapter 
VI of this report. 

The terminal listings are presented first, followed by 
the computer offline output. 
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repeat login nuc@pn0 

login 1621p44 400k 
ENTER PASSWORD: 

ENTER 4-DIGIT PROJECT NUMBER FOLLOWED BY 4-CHARACTER COST CENTER CODE 
0444C172 

READY AT 15.12.34 ON 09/16/78 
CMS Version 5.25 

CP 

link 2069p 191 192 
ENTER PASSWORD; 

SET TO READ ONLY 

ipl cms 

CMS Version 3.25 
login 191 

■ 

login 192 t»p ^ * 

T (192) R/0 
R; 

timeser 

EACH 2 SECONDS EXECUTION TIME IS INDICATED BY *. 

YOU HAVE ENTERED THE TIME SERIES EDITOR. - 

PLEASE RESPOND TO 'feSCH' UOERY 'WITH AN INPUT AT THE TERMINAL. 

ENTER ONLY THE FIRST LSTTER FOR A WORD RESPONSE. 

ENTER NUMERICAL VALUES VIA FORTRAN FORMAT. 

TYPE INTEGER VALUES (RIGHT JUSTIFIED)' FOR NAMES STARTING 

WITH I THRO N. TYPE FLOATING VALUES. WITH DECIMAL FOR ALL OTHERS. 

DO YOU WANT A LIST OP THE OPTIONS? 

y 



OPTION DESCRIPTION 

GENERATE GENERATE ANY ARIMA TIME ^ERIES ' 

AUTO CALCULATE AUTOCORRELATIONS, PAOTOS, MEAN AND VARIANCE 

PLOT —PLOT A TIME SERIES 

ESTIMATE CALCULATE MAX LIKELIHOOD ESTIMATES OF ARMA PARAMETERS 

DIFF DIFFERENCE A TI^E SERIES • 

FORECAST FORECAST FUTURE VALUES, '’CONSTRUCT CONFIDENCE INTERVALS 

TRANS --TRANSFORMS VALUES OF A TIME SERIES 

BOOTS DETERMINES ROOTS OF ARIMA CHARACTERISTIC EQUATION 

ZFORMAT ALTER DATA FILE TO FORMAT &Fl5;-6 

CMSWORK PERFORM CP/CMS COMMANTDS IN TIMESER EXEC 

SIMULATE SIMULATE NONSEASONAL TIME SERIES 

YESTSEAS CALCULATE INITIAL SEASONAL PARAMETERS 

WMARQRDT MARQUARDT SOLUTION FOR PARAMETER ESTIMATES 

XSUMSQ CALCULATE SON OP SQUARES FOR ARBITRARY PARAMETERS 



WOULD YOU LIKE MORE INFO? 

y . ■ 

ENTER OPTION YOU WANT INFO ABOUT. 

y 



YESTSEAS THIS PROGRAM CALCULATES INITIAL PARAMETER ESTIMATES 

FOR SEASONAL AND NONSEASONAL ARIMA MODELS, TO BE USED AS INPUTS TO 
THE WMARQRDT PROGRAM. 



DO YOU WANT INFO ABOUT ANOTHER OPTION? 

y 

ENTER OPTION YOU WANT INFO ABOUT, 

w 
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W-IARQRDT THIS PROGRAM CALCULATES NON-LINEAR LEAST SQUARES 

ESTIMATES OF BOX-JENKINS PARAMETERS FOR A GENERAL {SEASONAL 
OR NONSEASONAL) ARIMA model. it REQUIRES INITIAL PARAMETER 
ESTIMATES AS STARTING VALUES; THESE MAY BE CALCULATED USING 
PROGRAM VESTS EAS. 



DO YOU WANT INFO ABOUT ANOTHER OPTION? 

y 

ENTER OPTION YOU WANT INFO ABOUT. 

X 

XSUMSQ THIS PROGRAM CALCULATES THE SUM OF SQUARED RESIDUALS 

FOR ANY SET OF ARIMA (SEASONAL OR NONSEASONAL) PARAMETERS. 



DO YOU WANT INFO ABOUT ANOTHER OPTION? 
n 

*DO YOU WANT TO TRY A SESSION? 

y 



ENTER LETTER FOR OPTION YOU WANT, 
c 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 

WHEN FINISHED TYPE; 6GOTO -QUES 
stat 

P (191); 26 FILES; 235 REC IN USE, 61 LEFT (of 296), 79% FULL (2 CYL) 
cp q f 

FILES;- NO RDR, NO PRT, NO PUN 
listf * ft02f001 

FILENAME FILETYPE MODE NO. REC. DATE 
SERG FT02F001 PI 3 9/16 

SERC FT02F001 PI -5 • 9/16 

alter secg ft02f001 pi file ft02f001 pi 
sgoto -ques 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT. 

P 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

ENTER TITLE FOR PLOT 
series g data / undifferenced 
TIME SERIES PLOTS HAVE BEEN PRINTED OFFLINE 
DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT, 
a 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

AUTOCORRELATIONS 



0.948 


0.876 


0.807 


0.753 


0.714 


0.632 


0.663 


0.656 


0.671 


0.703 


0.743 

0.457 


0.760 

0.482 


0.713 

0.517 


0.646 

0.532 


0.586 

0.494 


0.538 


0.500 


0.469 


0.450 


0.442 



PARTIAL AUTOCORRELATIONS 
0.948 -0.229 0.038 0.094 0.074 


0.008 


0.126 


0.090 


0.232 


0.166 


0.171 -0.135 -0.540 -0.027 0.091 


0.025 


0.032 


0.07: 


0.048 


-0.046 



0.046 -0.100 0.052 0.048 -0.163 
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MEAN= 280.299 



VARIANCE * 14292.0 



ENTER TITLE FOR PLOTS. 

autos and pautos of series g data / undiffecenced 
YOOR AUTO AND PAOTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

DO YOU WANT TO GO AGAIN? 

y . 

ENTER LETTER FOR OPTION YOU WANT, 
t 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

DO YOU WANT/TO TRANSLATE THE ORIGIN; W«Z-SHIFT? 
n 

DO YOU WANT TO RESCALE THE VALUES: W»Z*SCALE? 
n 

DO YOU WANT A LOG TRANSFORMATION? 

y 

TRANSFORMATION IS W ( I ) »LOG(SCALE* (Z ( I )-SHIFT )-FACTOR) WHERE: 

SCALE* 1.00000 SHIFT*0.0 FACTOR-0.0 

YOUR ORIGINAL SERIES IS IN DATA FT02F001 Pi. 

YOUR TRANSFORMED SERIES IS IN FILE FT02F001 PI. 

DO YOU WANT TO PLOT THE TRANSFORMED VALUES? 

y 

EXECUTION BEGINS... 

* ENTER TITLE FOR PLOT 

plot of series g data / natural log transform 
TIME SERIES PLOTS HAVE BEEN PRINTED OFFLINS 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT, 
a 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

AUTOCORRELATIONS 



0.954 


0.899 


0.851 


0.808 


0.779 


0.756 


0.738 


0.727 


0.734 


0.744 


0.758 

0.498 


0.762 

0.506 


0.717 

0.517 


0.663 

0.520 


0.618 

0.484 


0.576 


0.544 


0.519 


0.501 


0.490 



PARTIAL AUTOCORRELATIONS 












0.954 -0.118 


0.054 


0.024 


0.116 


0.044 


0.038 


0.100 


0.204 


0.064 


0.106 -0.042 


-0.485 


-0.034 


0.042 


-0.044 


0.028 


0.037 


0.042 


0.014 


0.073 -0.033 


0.061 


0.031 


-0.194 













MEAN- 5.54218 VARIANCE - 0.193531 

ENTER TITLE FOR PLOTS. 

autos and pautos of series g data / natural log transform 
YOOR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

DO YOU WANT TO GO AGAIN? 

y 
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enter letter for option you want. 
d 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

IS YOUR TIME SERIES SEASONAL? 

y 

ENTER ORDER’ OF SEASONAL DIFFERENCING. 

1 

ENTER LENGTH OF SEASONAL PERIOD VIA 12. 

12 

ENTER NUMBER OF NONSEASONAL DIFFERENCES. 

1 

DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 

y 

EXECUTION BEGINS... 

AUTOCORRELATIONS 

-0.341 0.105 -0.202 0.021 0.056 0.031 -0.056 -0.001 0.176 -0.076 
0.064 -0.387 0.152 -0.058 0.150 -0.139 0.071 0.016 -0.011 -0.117 
0.039 -0.091 0.223 -0.018 -0.100 



PARTIAL AUTOCORRELATIONS 



-0.341 -0.013 
0.047 -0.339 
0.132 -0.072 



-0.193 -0.125 
-0.109 -0.077 
0.143 -0.067 



0.033 0.035 
-0.022 -0.140 
-0.103 



-0.060 -0.020 
0.026 0.115 



0.226 0.043 

0.013 -0.167 



MEAN«0.290920E-03 VARIANCE » 0.208604E-02 
ENTER TITLE FOR PLOTS. 

autos and pautos for series g s@data / 1 ns di££, 1 seas di££. In xform 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT. 

C 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 

WHEN FINISHED TYPE: &GOTO -QUES 
alter file £t03£001 pi save £t03£001 pi 
alter data £t02£001 pi serg £t02£001 pi 
stat 

P (191): 29 FILES; 267 REC IN USE, 29 LEFT (of 296), 90% FULL (2 CYL) 
erase file ftOBfOOl 
cp q f 

FILES:- NO RDR, NO PRT, NO PON 
igoto -ques 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT. 

y 

IS YOUR DATA IN FILE PT02F001? 

y 

EXECUTION BEGINS... 

ENTER LENGTH OP SEASON VIA FORMAT 12. 

12 

ENTER NUMBER OF NON-SEASONAL DIFFERENCES. 

1 

ENTER NUMBER OF SEASONAL DIFFERENCES. 

1 

ENTER NUMBER OP NON-SEASONAL AR PARAMETERS. 

0 
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ENTER NUMBER OF SEASONAL AR PARAMETERS. 

0 

enter number of non-seasonal ma parameters. 

^ENTER NUMBER OF SEASONAL MA PARAMETERS. 

1 

INITIAL PARAMETER ESTIMATES FOR MARQRT : 



THETA(l) - 0.390425 

THETAS(l) » 0.534397 

00 YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT, 
w 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

ENTER NUMBER OF NON-SEASONAL DIFFERENCES. 

1 

ENTER NUMBER OF SEASONAL DIFFERENCES. 

1 

ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 

0 

ENTER NUMBER OF SEASONAL AR PARAMETERS. 

0 

ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 

1 

ENTER NUMBER OF SEASONAL MA PARAMETERS. 

1 

ENTER LENGTH OF SEASON VIA FORMAT 12. 

12 

NOW INPUT YOUR INITIAL PARAMETER ESTIMATES r AS REQUESTED. 

ENTER NON-SEASONAL MA PARAMETER THETA(l). 

.390425 

ENTER SEASONAL MA PARAMETER THETAS (1). 

.534397 



CONVERGENCE HAS BEEN REACHED IN MAX{ 7, 7) ITERATIONS. 



SELECTED OUTPUT FOLLOWS: 



PARAMETER ESTIMATE STANDARD ERROR 



THETA(l) - 0.377152 0.819933D-01 

THETAS(l) » 0.572387 0.780252D-01 



MOVING AVERAGE CONSTANT: THETAO » 0.000291 



CHI-SQUARE STATISTIC FOR RESIDUAL LACK OF FIT » 29.713571 
DEGREES OF FREEDOM = 38 

PROBABILITY OF EXCEEDING STATISTIC * 0.829438 
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DO YOU WANT TO PLOT AUTO AND PAUTO OF RESIDUALS? 

y 

EXECUTION BEGINS... 

AUTOCORRELATIONS 



0.009 0.026 

0.026 -0.020 
-0.027 -0.014 



-0.129 -0.105 
0.013 0.035 

0.213 0.009 



0.078 0.077 

0.067 -0.130 
-0.042 



-0.036 -0.033 
0.052 0.015 



0.103 -0.050 
-0.093 -0.090 



PARTIAL AUTOCORRELATIONS 



0.009 0.026 

-0.022 0.018 
0.042 -0.014 



-0.130 -0.104 
0.046 0.001 

0.143 -0.025 



0.088 0.068 
0.053 -0.114 
0.013 



-0.072 -0.029 
0.069 0.029 



0.151 -0.057 
-0.129 -0.134 



MEAN=0.236005E-02 VARIANCE » 0.139713E-02 
ENTER TITLE FOR PLOTS. 

autos and pautos of residuals / (0,l,l)x(0, 1^1)12 model of series g 

YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

★ 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT. 

y 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

ENTER LENGTH OF SEASON VIA FORMAT 12. 

12 

ENTER NUMBER OF NON-SEASONAL DIFFERENCES. 

1 

ENTER NUMBER OF SEASONAL DIFFERENCES. 

1 

ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 

1 

ENTER NUMBER OF SEASONAL AR PARAMETERS. 

0 

ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 

1 

ENTER NUMBER OF SEASONAL MA PARAMETERS. 

1 

•k -k It 



INITIAL PARAMETER ESTIMATES FOR MARQRT ; 



PHI(l) =» 0.112699 

THETA(l) - 0.490883 

THETAS(l) - 0.533679 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT, 
w 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

ENTER NUMBER OF NON-SEASONAL DIFFERENCES. 

1 

ENTER NUMBER OF SEASONAL DIFFERENCES. 

1 

ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 

1 
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ENTER NUMBER OF SEASONAL AR PARAMETERS. 

0 

ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 

1 

ENTER NUMBER OF SEASONAL MA PARAMETERS. 

1 

ENTER LENGTH OF SEASON VIA FORMAT 12. 

12 

NOW INPUT YOUR INITIAL PARAMETER ESTIMATES, AS REQUESTED. 

ENTER NON-SEASONAL AR PARAMETER PHI(l). 

0.112699 

ENTER NON-SEASONAL MA PARAMETER THETA(l). 

0.490883 

ENTER SEASONAL MA PARAMETER THETAS (1). 

0.533679 

*** 

CONVERGENCE HAS BEEN REACHED IN MAX(10,10) ITERATIONS. 



SELECTED OUTPUT FOLLOWS: 



PARAMETER ESTIMATE STANDARD ERROR 



PHI(l) » 0.146667 0.224464D 00 
THETA(l) » 0.509816 0.195934D 00 
THETAS(l) » 0.573388 0.784896D-01 



MOVING AVERAGE CONSTANT: THETAO » 0.000248 



CHI-SQUARE STATISTIC FOR RESIDUAL LACK OF FIT » 28.302543 
DEGREES OF FREEDOM * 37 

PROBABILITY OF EXCEEDING STATISTIC » 0.847102 

DO YOU WANT TO PLOT AUTO AND PAUTO OF RESIDUALS? 

y 

♦EXECUTION BEGINS... 

AUTOCORRELATIONS 



-0.013 0.073 

0.033 -0.018 
-0.020 -0.023 



-0.100 -0.088 
0.014 0.027 

0.208 0.005 



0.077 0.073 

0.070 -0.133 
-0.036 



-0.031 -0.033 
0.051 0.008 



0.105 -0.054 
-0.092 -0.084 



PARTIAL AUTOCORRELATIONS 



-0.013 0.073 

-0.027 0.017 

0.051 -0.006 



-0.098 -0.096 
0.047 -0.004 
0.148 -0.019 



0.092 0.081 
0.053 -0.118 
0.003 



-0.064 -0.040 
0.055 0.031 



0.152 -0.050 
-0.123 -0.135 



MEAN=0.253052E-02 VARIANCE » 0.139373E-02 
ENTER TITLE FOR PLOTS. 

autos and pautos of residuals / (l,l,l)x(0,l, 1)12 model of series g 

YOUR AUTO AND PAUTO PLOTS ‘ HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

DO YOU WANT TO GO AGAIN? 

y 
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Ei^TER LETTER FOR OPTION YOU WANT, 
f ^ 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS.., 

IS YOUR SERIES SEASONAL? 

y 

ENTER SEASONAL DATA AS REQUESTED: 

• ENTER NUMBER OF SEASONAL DIFFERENCES, 

1 

ENTER NUMBER OF SEASONAL AR PARAMETERS, 

0 

ENTER NUMBER OF SEASONAL MA PARAMETERS. 

1 

ENTER LENGTH OF SEASON VIA FORMAT 12. 

12 

ENTER NONSEASONAL DATA AS REQUESTED; 

ENTER NUMBER OF NON-SEASONAL DIFFERENCES, 

1 

ENTER NUMBER OF NON-SEASONAL AR PARAMETERS, 

0 

ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 

1 

NOW INPUT YOUR PARAMETER ESTIMATES: 



ENTER MA CONSTANT TERM, THETAO. 

.000291 

ENTER NON-SEASONAL MA PARAMETER THETA (1). 

.377152 

ENTER SEASONAL MA PARAMETER THETAS (1). 

.572387 

ENTER MAXIMUM FORECAST LEAD TIME VIA FORMAT 12, 

30 

ENTER INDEX FOR PLOT ORIGIN VIA FORMAT 13. 

100 

ENTER INDEX FOR FORECAST ORIGIN VIA FORMAT 13, 

131 

ENTER SIGNIFICANCE LEVEL FOR CONFIDENCE INTERVALS. 

0.10 

WAS YOUR DATA TRANSFORMED IN THE TRANS PROGRAM? 

y 



DO YOU WANT BASIC OUTPUT AT THE TERMINAL? 

y 



THE LAST 10 WVEC VALUES: 

0.418999E 03 0.461001E 03 

0.606000E 03 0.507999E 03 

THE 30 FORECAST VALUES; 



0.398539^ 03 
0,469710E 03 
0.460421E 03 
0,520593E 03 
0,705666E 03 
0,526812E 03 



0.415184E 03 
0.543456E 03 
0.405213E 03 
0.507548E 03 
0.590927E 03 
0.502002E 03 



0.472000E 03 
0.461001E 03 



0.395402E 03 
0.618194E 03 
0.448016E 03 
0.528791E 03 
0.519089E 03 
0.587782E 03 



0.535001E 03 
0.390001E 03 



0.462697E 03 
0.626276E 03 
0.466864E 03 
0.611991E 03 
0,456978fi 03 
0.573220E 03 



0.622000E 03 
0.432002E 03 



0.450971E 03 
0.524292E 03 
0.444748E ?3 
0.696357E 03 
0.505397E 03 
0.597386E 03 
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THE 30 UPPER FORECAST CONFIDENCE LIMITS; 



0.399598E 


03 


0.416254E 


03 


0,396481E 


03 


0,463785E 


03 


0.452068E 


03 


0.470813E 


03 


0.544567E 


03 


0,619311E 


03 


0,627399E 


03 


0.525422E 


03 


0.461557E 


03 


0.406354E 


03 


0.449173E 


03 


0,468030E 


03 


0.445925E 


03 


0.521779E 


03 


0.508743E 


03 


0,529995E 


03 


0,613203E 


03 


0.697577E 


03 


0,706894E 


03 


0.592163E 


03 


0.520332E 


03 


0.458229E 


03 


0.506664E 


03 


0,528090E 


03 


0,503292E 


03 


0,589083E 


03 


0,574532E 


03 


0,598708E 


0^ 



ALPHA FOR THE CONFIDENCE LIMITS IS; 0.100 

* 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT, 

f 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

IS YOUR SERIES SEASONAL? 

y 

ENTER SEASONAL DATA AS REQUESTED; 

ENTER NUMBER OF SEASONAL DIFFERENCES. 

1 

ENTER NUMBER OF SEASONAL AR PARAMETERS. 

0 

ENTER NUMBER OF SEASONAL MA PARAMETERS. 

1 

ENTER LENGTH OF SEASON VIA FORMAT 12. 

12 

ENTER NONSEASONAL DATA AS REQUESTED; 

ENTER NUMBER OF NON-SEASONAL DIFFERENCES. 

1 

ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 

1 

ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 

1 



NOW INPUT YOUR' PARAMETER ESTIMATES; 



ENTER MA CONSTANT TERM, THETAO. 

.000248 

ENTER NON-SEASONAL AR PARAMETER PHI(l). 

0.146667 

ENTER NON-SEASONAL MA PARAMETER THETA (1). 

0.509816 

ENTER SEASONAL MA PARAMETER THETAS ( 1 ) . 

0.573388 

ENTER MAXIMUM FORECAST LEAD TIME VIA FORMAT 12. 

30 

ENTER INDEX FOR PLOT ORIGIN VIA FORMAT 13. 

100 

ENTER INDEX FOR FORECAST ORIGIN VIA FORMAT 13. 

131 

ENTER SIGNIFICANCE LEVEL FOR CONFIDENCE INTERVALS. 

0.10 

WAS YOUR DATA TRANSFORMED IN THE TRANS PROGRAM? 

y 



DO YOU WANT. BASIC OUTPUT AT THE TERMINAL? 

y 



158 



THE LAST 10 WVEC VALUES; 



0.418999E 


03 


0.461001E 


03 


0.472000E 


03 


0.535001E 


03 


0.622000E 


0. 606000E 


03 


0.507999E 


03 


0.461001E 


03 


0.390001E 


03 


0.432002E 


THE 30 FORECAST VALUES: 














0.398650E, 


03 


0.414875E 


03 


0.395125E 


03 


0 . 462366E 


03 


0.450651E 


0.469348E 


03 


0.543045E 


03 


0.617700E 


03 


0.625753E 


03 


0.523897E 


0.460071E 


03 


0.404899E 


03 


0.447739E 


03 


0.466360E 


03 


0.444325E 


0.520096B 


03 


0.507066E 


03 


0.528257E 


03 


0.611381E 


03 


0.695632E 


0.704905E 


03 


0.590336E 


03 


0.518566E 


03 


0.456512E 


03 


0.504960E 


0.526113E 


03 


0.501400E 


03 


0.587074S 


03 


0.572532E 


03 


0.596632E 


THE 30 UPPER 


FORECAST CONFIDENCE LIMITS: 








0.399709E 


03 


0.415945E 


03 


0.396204E 


03 


0.463452E 


03 


0.451744E 


0.470448E 


03 


0.544151E 


03 


0.618812E 


03 


0.626870E 


03 


0.525020E 


0.461199E 


03 


0.406032E 


03 


0.448887E 


03 


0.467517E 


03 


0.445491E 


0.521270E 


03 


0.508249E 


03 


0.529448E 


03 


0.612579E 


03 


0.696837E 


0.706117E 


03 


0.591556E 


03 


0.519792E 


03 


0.457745E 


03 


0.506207E 


0.527371E 


03 


0.502668E 


03 


0.588353E 


03 


0.573821E 


03 


0.597930E 


ALPHA FOR 


THE 


CONFIDENCE LIMITS IS; 0. 


100 









DO YOU WANT TO GO AGAIN? 

y 



ENTER LETTER FOR OPTION YOU WANT, 
f 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

IS YOUR SERIES SEASONAL? 

y 

ENTER SEASONAL DATA AS REQUESTED; 

ENTER NUMBER OF SEASONAL DIFFERENCES. 

1 

ENTER NUMBER OF SEASONAL AR PARAMETERS. 

0 

ENTER NUMBER OF SEASONAL MA PARAMETERS. 

1 

ENTER LENGTH OF SEASON VIA FORMAT 12. 

12 

ENTER NONSEASONAL DATA AS REQUESTED; 

ENTER NUMBER OF NON-SEASONAL DIFFERENCES. 

1 

ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 

0 

ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 

1 

NOW INPUT YOUR PARAMETER ESTIMATES; 



ENTER MA CONSTANT TERM, THETAO. 

.000291 

ENTER NON-SEASONAL MA PARAMETER THETA(l). 
.377152 

ENTER SEASONAL MA PARAMETER THETAS(l). 

,572387 

ENTER MAXIMUM FORECAST LEAD TIME VIA FORMAT 12. 
40 



03 

03 

03 

03 

o: 

03 

03 

03 

03 

03 

03 

03 

03 

03 
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ENTER INDEX FOR PLOT ORIGIN VIA FORMAT 13. 

075 

ENTER INDEX FOR FORECAST ORIGIN VIA FORMAT 13. 

120 

ENTER SIGNIFICANCE LEVEL FOR CONFIDENCE INTERVALS. 

0.10 

WAS YOUR DATA TRANSFORMED IN THE TRANS PROGRAM? 

y 



DO YOU WANT BASIC OUTPUT AT THE • TERMINAL? 

n 

* 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT, 
d 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 

IS YOUR TIME SERIES SEASONAL? 

y 

^ ENTER ORDER OF SEASONAL DIFFERENCING. 

ENTER LENGTH OF SEASONAL PERIOD VIA 12. 

12 

ENTER NUMBER OF NONSEASONAL DIFFERENCES. 

1 

DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 
n 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT. 

C 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 

VtfHEN FINISHED TYPE; &GOTO -QUES 
alter file ft02f001 pi Inserg ft02f001 pi 
alter file ft03f001 pi file ft02f001 pi 
&goto -ques 

DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT. 

X 

IS YOUR DATA IN FILE FT02F001? 

y 

EXECUTION BEGINS... 



IS YOUR SERIES SEASONAL? 

y 

ENTER LENGTH OF SEASON VIA FORMAT 12. 

12 

ENTER NUMBER OF SEASONAL AR PARAMETERS. 

0 

ENTER NUMBER OF SEASONAL MA PARAMETERS. 

1 

ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 

0 

ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 

1 



NOW INPUT YOUR INITIAL PARAMETER ESTIMATES, AS REQUESTED. 

ENTER NON-SEASONAL MA PARAMETER THETA(l). 

.377152 

ENTER SEASONAL MA PARAMETER THETAS (1). 

.572387 
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FOR THE GIVEN INPUT PARAMETERS AND MODEL 

THE RESIDUAL SUM OF SQUARES IS: 0.18192935D 00 



DO TOO WANT TO TEST DIFFERENT PARAMETER VALUES? 
n 

DO TOO WANT TO TEST A DIFFERENT MODEL? 

y 

IS TOUR SERIES SEASONAL? 

y 

ENTER LENGTH OF SEASON VIA FORMAT 12. 

12 

ENTER NUMBER OF SEASONAL AR PARAMETERS. 

0 

ENTER NUMBER OF SEASONAL MA PARAMETERS. 

1 

ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 

1 

ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 

1 

NOW INPUT TOUR INITIAL PARAMETER ESTIMATES, AS REQUESTED. 

ENTER NON-SEASONAL AR PARAMETER PHI(l). 

.146667 

ENTER NON-SEASONAL MA PARAMETER- THETA ( 1 ) . 

.509816 

ENTER SEASONAL MA PARAMETER THETAS ( 1 ) . 

.573388 



FOR THE GIVEN INPUT PARAMETERS AND MODEL 

THE RESIDUAL SUM OF SQUARES IS: 0.181643730 00 



DO YOU WANT TO TEST DIFFERENT PARAMETER VALUES? 
n 



DO YOU WANT TO TEST A DIFFERENT MODEL? 
n 

DO YOU WANT TO GO AGAIN? 
n 

CONTROL RETURNED TO CMS 

R; 
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APPENDIX D 



PROGRAM LISTINGS 

This Appendix contains program listings of all programs 
resident in the Time Series Editor. It does not include 
listings of subroutines used from other libraries, such 
as IMSL routines. 
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APPENDIX E 



LISTING OF DATA SETS RESIDING 
IN THE TIME SERIES EDITOR 

This appendix contains listings of the data sets now 
residing in the Time Series Editor. The data consists of 
Box and Jenkins' [Ref. 4] time series C [Ref. 4, p. 528] 
and time series G [Ref. 4 , p. 531], as well as a data set 
containing monthly Monterey, California, rainfall data 
from January 1931 through December 1976 . 
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